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Synopsis 


Hydrogenated amorphous silicon(a-Si:H) is the most widely studied material amor 
disordered semiconductors. a-Si:H thin fihns made by the glow discharge of silane(I 
have a low (10^® — density of localised states(DOS) at the Fermi ! 

This allows the material to be doped n-type or p-type at will. This is usually dor 
mixing SiH 4 with phosphine(for n-type) or diborane(for p-type) during the prepari 
by glow discharge. a-Si:H can also be doped interstitially by in diffusion of lith 
Thus, the material finds use in many devices like photovoltaic cells, photoconduci 
thin film transistors etc. 

DOS in the mobility gap of the bulk of a-Si:H is the most significant param 
that decides its suitability for device applications. Also, the density of surface state 
important. Various experiments like field effect, thermally stimulated currents, ca 
citance, deep level transient spectroscopy, space charge limited currents and sub | 
absorption etc have been used for the determination of DOS in a-Si:H. Each metl 
has its advantages and limitations. Optical absorption in the sub gap region[a(/ii/)] ] 
been used for the determination of bulk states with the advantage that the informati 
about the shallow as well as the deep states can be obtained. 

Two techniques, namely, constant photocurrent mecLSurements(CPM) and phol 
thermal deflection spectroscopy (PDS) are used to measure the sub gap absorption 
a-Si:H. One generally finds, however, that the two techniques yield different values 
a(/n'), on the same sample. This is expected, since CPM is a photocurrent measur 
ment and is not sensitive to the surface states. In a-Si:H, the surface states are larj 
in number and the carriers generated in the surface region recombine quickly wither 
contributing to the photocurrent. In addition, some of the bulk states may also nc 
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contribute to CPM if the transitions they participate in, have localised states as the 
final states; On the other hand, PDS is a thermal effect and is sensitive to all the 
transitions involving the bulk as well as the surface states. 

In the literature, the density of surface states has been obtained by measuring oc[hi/) 
on identical a-Si:H films of varying thicknesses, using PDS. However, making such 
samples is practically difficult and gives large uncertainty in the measured density, pf 
surface states. 

Jn the present work, a measured by CPM and PDS on a given sample of given 
thickness is analysed to obtain the distribution and density of bulk states in the entire 
mobility gap as well as surface states. This method obviates the difficulties associated 
with' the 'preparation of identical samples of varying thicknesses. In addition, the dual 
beam constant photocurrent measurements(DBCPM), are used to obtain the correlation 
energy between electrons. Further, using DBCPM we also obtain the shape of the mid- 
gap states, which is not possible using CPM or PDS. We have applied this technique of 
measuring the bulk and surface states in undoped a-Si:H and a-Si:H doped with various 
amounts of phosphorous (P). Our results are in qualitative agreement with the literature. 
We also report the measurements of a{hi/) on lithium(Li) doped samples. This is 
interesting because Li is interstitial dopant in a-Si:H, unlike P which is substitutional. 
Since, the shape of the surface states in a-Si:H is not known, various shapes of surface 
states like a delta function, a rectangular and a Gaussian are tried. We find that the 
calculated density of surface states is not very sensitive to the shape chosen. Further, 
the density of bulk and surface states in P and Li doped a-Si:H are found to be similar. 

For device fabrication, a-Si:H may be subjected to various surface treatments like 
oxidation, etching with hydrofluoric acid and exposure to various plasma. These may 
change the density of states in the bulk and at the surface. These changes can give 
useful information for device fabrication and have not been studied in detail, so far. We 
have measured changes in the density of bulk and surface states after various surface 
treatments, and light soaking using sub-gap absorption measurements. We find that 
light soaking(Staebler Wronski effect) is mainly a bulk effect and does not change the 
density of surface states much. Native oxide etching, using hydrofluoric acid, is found 
to increase the density of surface states, whereas, treatment with hydrogen plasma at 
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room temperature decreases the density of surface states. This might explain why 
treating a-Si:H with hydrogen plasma, in device fabrication, yields better results. 

The potential fluctuations play an important role in determining the electronic prop- 
erties of a-Si:H. Several experiments like NMR, show that hydrogen is distributed non- 
imiformly in a-Si:H. A non-uniform distribution of hydrogen is likely to give rise to 
potential fluctuations in a-Si:H. In addition, presence of charge centers, charged and 
neutral dangling bonds and the dopants, if present, may also contribute to them. Over- 
hof and Beyer have used transport measurements(conductivity and thermopower) to 
obtain the width of these potential fluctuations in a-Si:H. We have estimated the width 
of potential fluctuations using sub-gap absorption measurements, on P doped a-Si:H. 
We find that the width of potential fluctuations increases with the increasing doping 
concentrations. In addition, we have also measured the width of long range poten- 
tial fluctuations by transport measurements on the same samples, using the analysis 
of Overhof and Beyer. The width of the potential fluctuations, estimated using op- 
tical absorption measurements, is foimd to be about a factor of 20 more than that of 
the long range potential fluctuations, estimated using transport measurements. These 
results are explained qualitatively, by arguing that in transport measurements, carri- 
ers are likely to overlook the sharp fluctuations in the potential because of quantum 
mechanical tunneling of carriers, as these fluctuations are short ranged. On the other 
hand, sub gap absorption measurements are local processes and are likely to be sensit- 
ive to the short range potential fluctuations, which may have a larger width. We also 
explain qualitatively the difference in a measured by CPM and PDS using potential 
fluctuations. 

Organisation of the thesis 

The present thesis has seven Chapters. Chapter 1 gives a brief introduction to the sub 
gap absorption and the potential fluctuaitions in a-Si:H and the motivation behind the 
present work. The sample preparation, doping and initiad characterisation of the films 
are given in Chapter 2. Details of the constant photocxirrent measurement (CPM) and 
photothermal deflection spectroscopy(PDS) set ups are adso contained in this Chapter. 
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In Chapter 3, we describe the determination of bulk density of states in the entire 
mobility gap using CPM and PDS measurements. The dual beam CPM results used 
to determine the shape of the mid gap states, are also discussed in this Chapter. 

a.{hu) measured by CPM and PDS have been used to determine the surface states in 
Chapter 4. Results on undoped a-Si:H and a-Si:H doped with varying phosphorous and 
lithium concentrations are also presented in this Chapter. Chapter 5 gives the results 
on the determination of density of surface states after various surface treatments. In 
Chapter 6, we explain the difference in a(hi/) as measured by CPM and PDS, in terms 
of the potential fluctuations present in a-Si:H. The width of the potential fluctuations 
estimated by transport measurements is compared with those estimated by the optical 
absorption measurements. Finally, Chapter 7 summarises our results and lists the 
major conclusions of the thesis. Some unanswered questions and directions for further 
research are suggested. 
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Chapter 1 


Introduction 


Hydrogenated amorphous silicon (a-Si:H) is the most promising tetrahedrally bonded 
amorphous semiconductor, being used for many devices like photovoltaic cells [1], pho- 
toconductors, thin film transistors [2] and many more [3]. There are many proposed 
applications of the material [3] and therefore it is important to study and solve the 
persisting problems and issues. 

The difference between the amorphous semiconductors and their crystalline counter- 
part is the lack of long range periodicity of atomic structure in the former. But the short 
range periodicity is preserved even in the amorphous semiconductors [4]. This results 
in a similar electronic structure like the band gap, as in the crystalline semiconductors. 
Thus, amorphous semiconductors have bands consisting of extended states separated by 
a mobility gap which may have localised states coming from disorder [5]. This disorder 
is the result of variation in the bond angle and the bond lengths. This results in the 
taihng of the band edges giving rise to the localised states in the tail region. This is 
known as the Anderson locahsation [5]. In addition, amorphous semiconductors may 
also have coordination defects (like dangling bonds) resulting in some more localised 
states in the middle of the mobihty gap. The density of localised states in the gap is 
usually large (?s 10^° [6] in an amorphous semiconductor. 

Because of the high density of locahsed states (DOS), near the Fermi level (E;?). 
the effect of doping in most amorphous semiconductors can not be observed, as their 
Fermi level is pinned. Spear and LeComber [7] discovered in 1972 that unlike other 
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amorphous semiconductors, a-Si:H prepared by the glow discharge of silane has a low 
density of localised states at E^. Also, the dark conductivity can be changed by several 
orders of magnitude by gas phase doping either n type or p type [8,9]. It is believed that 
hydrogen in a-Si:H is crucial in achieving low density of states at Ei?. Further studies 
revealed that hydrogen not only act as the bond terminators, reducing the density of 
dangling bond states, but also releases the strain in the network, reducing the bond 
length and bond angle fluctuations. 

1.1 Bulk and Surface States 

Presence of disorder and defects in a-Si:H give rise to the localised electronic states in 
the mobility gap. These states signiflcantly affect the properties of the material. Also, 
doping changes the density, nature and the energy levels of these states. Therefore, the 
information about the localised states is of great importance. 

The mobility gap in a-Si;H, contains a continuum of states. We get exponentially 
decreasing shallow states near the band edges, and deep states arising from the coordin- 
ation defects. Fig. 1.1 shows the schematic of the standard model of bulk localised states 
in undoped a-Si:H, as proposed by Street et al [10]. Doping changes the distribution 
and density of gap states [11]. 

It increases the bond distortions and also increases the deep gap states. The nature 
and the position of the deep gap states also change. On n-type doping, the singl}' 
occupied neutral dangling bond states (D°) change to doubly occupied states D“, and 
we get only one Gaussian peak for the deep states as compared to two peaks in the case 
of undoped a-Si;H [10]. .. 

Localised states in the band gap have been extensively studied by Luminescence 
[10]. electron spin resonance [12], field effect measurements [13], thermally stimulated 
currents [14], space charge limited currents [15], capacitance techniques [16] and sub 
gap absorption measurements [17-20] etc. Sub-gap absorption is a fundamental process 
that measures the density and the energ}^ position of all the states in the band gap, thus 
providing a convenient tool for DOS characterisation. However, certain assumptions are 
necessary to obtain information of DOS from the sub gap absorption data. The matrix 
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Figure 1.1: Bulk localised states consisting of exponential lirbacli edges and Gaussian 
mid-gap states. E„ and are band edges. Ej^ is the Fermi level. D° and D" are 
neutral and charged dangling bond states. U is the correlation energy [after reference 



elements of the optical transitions are not known. These are generally assumed to be 
independent of the photon energy, i.e. the matrix element for the transitions between 
the localised to the extended states are assumed to be the same as that from extended 
to extended states. Also, in order to get DOS from sub-gap absorption, the shape of 
the distribution of DOS has to be known, in advance, and has to be guessed. Various 
shapes of DOS have been suggested [10,21]. So far, the issue has not been resolved. 

The smallest absorption coefficient (a) that can be measured comfortably using the 
conventional transmission and reflection techniques is 50-100 cm-^ for a 1 ^m thick 
a-Si:H film. The a associated with the transitions involving the deep states in a-Si:H 
can be smaller by a factor of about 100. In an attempt to measure extremely small 
Q, two techniques, namely, constant photocurrent measurement (CPM) [17,18] and 
photothermal deflection spectroscopy (PDS) [19,20] are employed. We describe the 
CPM and PDS set ups in Chapter 2. 

a{hu) measured by the CPM and PDS on the same sample, are usually not equal, 
Opds being larger than a^pm- Fig. 1.2 illustrates this for an undoped a-Si:H sample 
(#18). The difference arises because CPM is not sensitive to the optical transitions 
involving a part of the bulk states and full surface states [22,23], whereas, PDS is 
sensitive to all the transitions. 

We analyse acpm.{hu) and apds{hu) data to get the bulk states in the entire mobility 
gap (Chapter 3), and surface states (Chapter 4), for undoped and doped a-Si:H. This 
method of measuring the density of surface states (N,), has the advantage, that it needs 
only one sample. The usual method for measuring N, needs a number of identical 
samples of varying thicknesses [22,24,25]. This is difficult to achieve, practically. 

In addition to jihosphorous (P) incorporation, a-Si:H can also be doped n-type by 
the in-diffusion of lithium (Li). However. Li is interstitial donor [27], whereas, P is sub- 
stitutional. Also, because Li doping is done at 100-150°C, the hydrogen concentration 
is likely to remain the same as in undoped a-Si:H. Therefore, one can study the effect 
of the dopant alone. It is interesting to compare the DOS in P doped 
a-Si:H [a-Si:H(P)] and Li doped a-Si:H [a-Si:H(Li)]. 

Also, It is worthwhile to see if some of the issues that remain unresolved for a-Si:H(P ) 
remain unanswered for a-Si:H(Li) also. One such controversy is about the position of 
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Figure 1.2. Sub-gap absorption coefficient measured as a function of photon energy (hi/) 
using CPM (solid circle), PDS (open circles) and transmission (crosses) [after reference 




D in a-Si;H(P). Street et al [10] suggest that the position of D" in undoped and in 
P doped a-Si:H are the same. This is called the ’’Basic Model”. On the other hand, 
Kocka et al [27], based on sub-gap absorption measurements, suggest that on doping, 
the dangling bonds form intimate pair with the dopants, resulting in the shift of the 
energy position of D“ away from the conduction band edge. Now, Li may or may not 
form the intimate pair with dangling bonds. Therefore, it is interesting to study the 
position of dangling bond states in a-Si:H(Li). 

In a-Si:H based device fabrication, a-Si:H has to undergo various surface treatments. 
It has been found that treatment with hydrogen plasma improves the quality of devices 
having a number of interfaces, like solar cells. A possible reason for this may be that the 
treatment decreases the number of surface states. However, in the literature the density 
of surface states have not been measured before and after these surface treatments. We 
report estimation of after light soaking and various surface treatments like etching, 
oxidation and exposure to various plasma etc (Chapter 5). 


1.2 Potential Fluctuations 

a-Si:H grown by glow discharge is heterogeneous [28-30], Structural investigations of 
the a-Si:H Hms show that the bonded hydrogen is distributed non-uniformly in clusters, 
even in the best quality material [31], The heterogeneities give nse to potential fluc- 
tuations. in addition, potential fluctuations in a-Si:H may also arise because of the 
variations m the bond angle and bond length, non-uniform distribution of the neutral 
and charged dangling bonds as well as the dopants (if present). It is difficult to handle 
the problem of potential fluctuations, in a fundamental way. Nevertheless, they are im- 
portant as they must influence the electronic properties of a-Si:H in a profound maoo e . 
[32j. One such property is the degradation of electronic properties of a-Si:H on light 
soaking (Staebler-Wronski effect) [33]. It is found that additional dangling bonds are 
created [34], on light soaking. Other photo induced structural changes have also been 

suggested [32]. Also, there are reports of increase in the width of potential fluctuations 

[35-37], on light soaking. 

In a-Si:H, the range of potential fluctuations from structural disorder, is likely to 
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be of the order of inter-atomic distancae [38] (« 3A). In transport measurements, 
carriers are required to move through the potential fluctuations. Therefore, carriers 
can’t see the sharp changes in the potential, because they can tunnel through them, 
quantum mechanically, therefore, the transport measurements are sensitive to the long 
range potential fluctuations [35,36). Fig 1.3 shows the in phase part of the long rangl 
potential fluctuations, associated with non uniform distribution of charged centers or 
region which may be related to impurity and structural disorder [28], Local density 
fluctuations or compositional variations are expected to produce spatial variations of 
the mobility gap and have not been shown in the figure, for clarity. It should be noted 
that the hnes drawn represent the local potential average over a region of the order of 
coherence length of the electron and not the true local potential. In such fluctuating 
mobility edges, spatial average of conductivity is dominated by highest energy regions 
of the potential fluctuations where the conductivity is lowest. On the other hand, the 
thermopower measures the average energy over the complete transport path and the 
low and high energy regions of the potential fluctuations contribute more equally. The 
difference (Eg) m the conductivity and thermopower activation energies is, therefore, a 
rough measure of the magnitude of potential fluctuations. 

Overhof and Beyer [35,36] have done model calculations and have related to the 
width of long range potential fluctuations to Eg quantitatively. 

Whereas, transport measurements are sensitive to the long range potential fluctu- 
ations, optical absorption measurements are likely to be sensitive to the short range 
potential fluctuations, of the order of interatomic distances. This is because optical 
absorption are vertical transitions in the real space (Fig. 1.3). Therefore, an analysis 
of sub-gap absorption measurements using potential fluctuations, is likely to give the 
width of short range potential fluctuations in contrast with the transport measurements, 
which give the information about the long range potential fluctuations. 

In a-Si:H, oc{hi^) has Urbach edge region (region II) and plateau region (region III. 
see Fig. 1.2). Absorption in the Urbach edge region is attributed to bond length and 
bond angle variations, whereas region II is attributed to the absorption involving the 
deep gap states. John et al [38] have done quantum mechanical calculations to obtain 
exponential Urbach edge starting from Gaussian potential fluctuations. They relate the 
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Figure 1.3. Long range potential fluctuations as might be present in a-Si:H(Ei/ and 
Ec represent the valance and conduction band percolation edges, respectively. Ep is 
the Fermi level. Note that short range potential wells are not shown. The part of the 
potential fluctuations which causes changes in E, is also not shown (after reference 28). 
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width of the potential fluctuations (Vr) to slope of Urbach edge (E„) [38], 

1.3 Motivation for the Present Work 

In the previous sections, we very briefly described the current status of this field and 
the problems and puzzles that need to be solved. We have undertaken a study of some 
of these issues and they are the main motivations of the present work. In the following 
we describe the salient features of the problems that have motivated this study. 

La the literature, aaprn.{hu) and aj,tis{hv) have been used to obtain the density of 
localised states only in a part of the mobility gap. In particular, the slope of the 
conduction band Urbach edge has not been obtained using sub-gap absorption meas- 
urements. Further, for the DOS determination, the shape of DOS has to be guessed, 
in advance. Several different shapes have been proposed. We have used dual beam 
CPM (DBCPM) data and tried to guess the shape of DOS (Chapter 3). In addition, 
the value of the correlation energy (U), is not known precisely, although it is agreed 

that U IS not negative [39]. We have attempted to measure U, using DBCPM (Section 
3.1.1.4). 

Li IS an interstitial dopant m a-Si:H. It is observed that on doping with phosphorous 
from the gas phase, less hydrogen gets incorporated in the material. Since, Li diffusion 
is done at lower temperatures (less than ISO^C), the hydrogen concentration in the 
material is expected to remain the same as that in undoped samples. This allows one 
to study the effect of doping in a-Si:H(Li) without worrying about changing hydrogen 
concentration. In particular, it is interesting to compare a-Si:H(P) and a-Si:H(Li) for 
their bulk as well as surface states. 

The existing methods of determining density of surface states using sub-gap absorp- 
tion technique, presents practical difficulties of preparing a number of identical films. 
Other methods need perfectly uniform film and can be used only for device qualit}^ 
materials. Thus a technique, free from the above mentioned difficulties, is desirable 
which can estimate bulk states as weU as the surface states. Our analysis, described in 
Chapter 3 and 4, obviates these difficulties. 

Potential fluctuations in a-Si:H are not studied much. For example, difference 
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between a{hi/) measured by CPM and PDS have not been explained using potential 
fluctuations. Also, it is observed that transport measurements see potential fluctuations, 
which is some kind of average of the short range potential fluctuations. Optical absorp- 
tion is sensitive to the short range potential fluctuations. Therefore, width of potential 
fluctuations measured by transport and optical measurements on the same samples may 
give how the potentials, as seen by transport measurements, are averaged. These and 
other issues have been taken up in Chapter 6. 

1.4 Organisation of the Thesis 

The next Chapter describes the preparation of undoped and doped a-Si:H films by glow 
discharge method and doping with phosphorus and lithium. The characterisation of 
these films for amorphous nature, dark conductivity as a function of temperature, room 
temperature photoconductivity and thermopower is discussed in brief. Experimental 
set-up for the transmission measurement for obtaining absorption coefficient (a) in the 
region 2.5eV> hz/ > 1.6eV, thickness and band gap of a-Si:H thin films is discussed. The 
details of the instrumentation for computer controlled measurements of sub-gap a using 
constant photocurrent measurement (CPM) and photothermal deflection spectroscopy 
(PDS) techniques are also given in this Chapter. 

Chapter 3 describes the principle and analysis of a^ds{hu) and acj,m{hu) data. The 
technique adopted for the determination of the bulk states in the entire mobility sap 
in undoped and n type doped a-Si:H is discussed. Further, the details of dual beam 
constant photocurrent measurement (DBCPM) to obtain the correlation energy U and 
the shape of the deep gap states are discussed in this Chapter. Results on undoped a- 
Si:H, a-Si:H doped with various quantities of Li and P are discussed. Also, the position 

of D in a-Si:H(Li) have been discussed. Finally we give summary and conclusions of 
this Chapter. 

In Chapter 4, we apply acpM and ap^s data for the measurement of densit\- of 
surface states. Again results on the samples, for which the bulk states are determined in 
Chapter 3, have been discussed. The effect of detailed distribution of the surface states 
have been discussed. Finally, we give the summary and conclusion of this Chapter. 
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Low density of surface states are crucial from the device point of view. In Chapter 
5, we report the results on the measurement of density of surface states after various 
surface treatments like light soaking, oxide etching, oxidation using concentrated H 25'04 
and annealing and treatment with room temperature nitrogen and Hydrogen plasma. 
These experiments also study the sensitivity of the technique for the measurement of 
density of surface states proposed in Chapter 4. 

Chapter 6 describes how width of short range potential fluctuations and the 

long range potential fluctuations (A) can be obtained using sub-gap absorption and 
transport measurements respectively. Sub-gap absorption and transport measurements 
are performed on the same samples and and A are calculated. An empirical 

relation between Vrms and A is found. An attempt have been made to explain this 
empirical relationship. Further, a more comprehensive distribution of the potential 
fluctuations have been proposed, which explains both sub-gap absorption data in the 
entire region measured by us. 

Chapter 7 summarises the results and the conclusions of the present study. The 
directions to carry out further work which helps in understanding the bulk and surface 
states and the potential fluctuations in a-Si:H are also given in the Chapter. 


11 



Bibliography 

1. S. Guha, J. Non. Cryst. Solids, 198 -200, 1076 (1996). 

2. S. Tomiyama, T. Ozawa, H. Ito and T. Nakamura, J. Non. Cryst. Solids, 198 
-200, 1087 (1996). 

3. R. A. Street, S. Nelson, L. Antonuk and V. Perez mendez, MRS Symp. 192 441 
(1990). 

4. N. F. Mott and E. A. Davis, Electronic processes in Non-Crystalline Materials, 
(Clarendon Press, Oxford, 1979), Chapters 2 and 6. 

-5. S. R. Elliot, Physics of Amorphous Materials, (Longman, UK, 1990), Chapter 3. 

6. H. Overhof and P. Thomas, Electronic Transport in Hydrogenated Amorphous 
Semiconductors, Volume 114, (Springer- Verlag, Berlin, 1989). 

7. W. E. Spear and P. C. LeComber, J. Non. Cryst. Solids, 8-10. 727 (1972). 

5. W. E. Spear and P. C. LeComber, Solid State Comm. 17, 1193 (1975). 

9. W. E. Spear and P. C. LeComber, Phil. Mag. 33, 935 (1976). 

10. R. A. Street, J. Zesch and M. J. Thompson, Appl. Phys. Lett. 43(7), 672 (1983). 

11. .J. Kocka, J.Non.Cryst. Solids, 90, 91 (1987). 

12. .J. D. Cohen, J. P. Habrison and K. W. Wecht, Phys. Rev. Lett. 48, 109 (1982). 

13. W. E. Spear and P. C. LeComber .J. Non. Cryst. Solids, 8-10. 727 (1972). 

D. S. Misra, Ph.D Thesis, Determination of the density of localised states in 
well characterised thin films of hydrogenated amorphous silicon prepared by glow 
discharge”, IIT Kanpur, (Unpublished) 1984. 

15. S. -A.shok and S. J. Fonash, IEEE. Electron Devices, ED-1, 200 (1980). 

16. P. Viktorovitch and D. Moddel, J. .4ppl. Phys. El, 4847 (1980). 


12 



17. M. Vanecek, J. Kocka, J. Stuchlik, Z. Koicek, 
Energy Materials 8, 411 (1983). 


0. Stika, and A. Triska, Solar 


18. A. K. Sinha. M. Malhotra, S. Kvunar. E. Bhatt^Earya and S. C. Agaxwal. Ind. 
J. Pure & Appl. Phys. 31, 548 (1993). 

19. W. B, Jackson, N. M. Amec. A. C. Boccara, and D. Fournier, AppI.Optics 20, 
1333 (1981). 


20. S. C. Agarwal, J. S. Payson and S. Guha, Phys. Rev. B 36, 9348 (1987); A. K. 
Sinha and S. C. Agaiwal, Ind. J. Pure k Appl. Phys. (1998) (in press). 

21. M. Hack and M. Shur, J. Appl. Phys. 54 , 5858 (1983); S. Guha, J.Non.Cryst. 
Solids, 77&:78, 1451 (1985). 

Z. E. Smith, V. Chu, K. Shepard, S. Aljishi, D. Slobodin, J. Kolodzey, S. Wagner, 
and T. L. Chu, Appl.Phys.Lett. 50, 1521 (1987). 

23. A. K. Sinha and S. C. Agarwal, Phil. Mag. B (1998), (In Press). 

24. H. Curtnis, N. Wyrsch and A. V. Shah, Electronics Letters 23, 228 (1987). 

25. M. Favre, H. Curtnis and A. V. Shah, J.Non-Cryst. Solids, 97 & 98, 731 (1987). 

26. W. Beyer and R. Fischer, Appl. Phys. Lett. 31, 850 (1977). 

27. J.Kocka, M. Vanecek and A.Triska. Amorphous silicon and Related Materials (ed. 
H. Fritzsche), Vol. lA, (World Scientific, Singapore, 1989), page 296. 

28. H. Fritzsche, J. Non-Cryst. Solids, 6, 49 (1971). 

29. S. C. Agarwal, Ind. J. Pure and .^ppl. Phys., 34, 597 (1996); S. C. Agarwal, 
Bull. Mater. Sci. 18, 669 (1995). 

30. E. N. Economou, C. M. Soukouhs, M. H. Cohen and A. D. Zdetsis, Phys. Rev. 

B 31, 6172 (1985); M. H. Cohen, in Random Media, IMA Vol. 7, (ed. G. 
Papanicolaou), Springer- Verlag xNew York (1987) page 89. 


13 



31. R. A. Street, Hydrogenated Amorphous Silicon (Cambridge University, 1991) 
Chapter 7. 

32. H. Fritzsche, Solid State Comm. 94, 953 (1995). 

33. D. L. Staebler and C. R. Wronski, Appl. Phys. Lett. 31, 32 (1977). 

. 34. H. Derch, J. Stuke and J. Beichler, Appl. Phys. Lett. 38, 456 (1981). 

35. H. Overhof and W. Beyer, Phil. Mag. B 43, 433 (1981). 

36. H. Overhof and W. Beyer, Phil. Mag. B 47, 377 (1983). 

37. P. Agarwal and S. C. Agarwal, J. Appl. Phys. (1997). 

38. S. John, C. M. Soukoulis, M. H. Cohen and E. N. Economou, Phys. Rev. Lett. 
57, 1777 (1986). 

39. See for example, J. Z. Liu, G. Lewen, J. P. Conde and P. Roca i Cabarrocas, J. 
Non-Cryst. Solids 164 Sz 165, -383 (1993). 



Chapter 2 


Experimental 


This Chapter describes the preparation of both undoped and doped hydrogenated 
amorphous silicon (a-Si:H) thin films and their initial characterisation for quality. Fur- 
ther, the details of the set ups for transmission, constant photocurrent measurement 
(CPM), photothermal deflection spectroscopy (PDS) and dual beam constant photo- 
current measurement (DBCPM) are also given. 

2.1 Preparation of a-*Si:H Films 

Hydrogenated amorphous Silicon (a-Si:H) is prepared by the decomposition of silane 
{Si Hi) gas in a mixture with argon or hydrogen in rf glow discharge apparatus [1]. 
Many variations of the technique have been tried to tailor the film properties [2,3]. Fig. 
2.1 shows the glow discharge reaction chamber made of stainless steel. A capacitor 
geometry is used. The substrates are loaded on the top electrode (anode), which is 
heated externally using a heater as shown in the figure. Prior to the deposition, the 
system is evacuated to the base pressure of fa 10~® torr using c. turbo molecular pump. 
The system is purged with hydrogen and substrates are cleaned using hydrogen plasma. 
Silane with the carrier gas is then introduced in the reaction chamber, and an electric 
field is applied between the cathode and the anode (Fig. 2.1). Silane molecules are 
dissociated through collisions with the fast electrons in the plasma, into neutral radicals, 
atoms and ions. These atoms and ions travel to the growing surface through a gas 
phase diffusion process, resulting in amorphous silicon thin films onto Corning 7059 
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Table 2.1: Optimum parameters for depositing undoped a-Si:H, using rf glow discharge 


Ga^ Composition 

Ratio 

SiH4 : H2 

10:90 

Flow Rate 

20 seem 

Power Density 

40mW/cm^ 

Deposition Rate 

1 A/sec. 

Substrate Temperature 

300°C 

Pressure of 


deposition 

1 torr 


glass substrates. More details of the deposition system are given elsewhere [4]. A multi 
chamber glow discharge system is described in ref. [5]. 

The quality of a-Si:H films depends upon the deposition conditions, like ratio of gases 
in the gas mixture, flow rate of gases, substrate temperature, pressure of deposition, 
absorbed power density etc. A combination of these parameters for obtaining good 
quality material varies from system to system. Table 2.1 gives the optimum deposition 
parameters for our system. Some of the films for this study have also been made outside 
optimum conditions to vary the density of defect states. 

2.1.1 Substrates 

Corning 7059 glass slides (25mmx25mmx0.5mm) are used as substrates for optical 
and electrical measurements. These are alkali free and are transparent in the region of 

interest (0.6eV to 4.0eV). For IR spectroscopy polished high resistivity silicon wafers 
are used. 

The 70o9 glass substrates are first washed in liquid detergent, rinsed in distilled 
water and then deionised water. These are then cleaned ultrasonically in acetone and 
finaUy degreased in vapours of isopropyl alcohol. Silicon substrates are degreased and 
then etched with 15% solution of HF. Substrates are then rinsed in deionised water and 
finally are exposed to the vapour of isopropyl alcohol. 
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Figure 2.1; Sdiematic diagram of the reaction chamber used for preparation of a-Si:H. 
Here, GV is the gate valve, B and PG are beratron and penning gauges. A, H and S 
are anode, substrate heater and substrate respectively. Cl and CL represent insulators 
and SS cyhnder. D is the view port, C the cathode and TC the thermocouple (after 
reference 5). 
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2.1.2 Doping 

For phosphorous (P) doping, SiH 4 is mixed with PH 3 in appropriate amounts. Phos- 
phorus acts as a substitutional donor. For interstitial doping [ 6 ], both vacuum evapor- 
ation of lithium (Li) and ion implantation of Li have been used [7,8], with success. 
have used vacuum evaporation of Li. 

Li pieces, covered with methanol in a tungsten boat, are loaded in a vacuum evap- 
oration system. Methanol is put to protect Li from air. a-Si.'H thin film is kept at an 
elevated temperature of 100 — 150°C and about 100 A thick Li is evaporated on top of 
it at a base pressure of about 10"® torr. The sample is then annealed at a temperature 
100 - 200 °C for 1-2 hours. Higher temperature and longer times of annealing give higher 
concentration of diffused Li in a-Si:H. The excess Li is removed from the surface by 
etching the film with 10 % HF. The sample is again annealed at about 150°C for 6 
hours to make the diffusion as uniform as possible. More details are given elsewhere 
[7.8]. 

2.1.3 Electrodes 

Xichrome and graphite (Aqua dag) are used as electrodes. It is found that both make 
ohmic contacts for electric fields of interest up to 600 V /cm. About l^m thick nichrome 
strips (ks 1 cm. in length and « 2 mm apart) are evaporated through a suitable mask 
on a-Si:H film deposited on 7059 glass substrate. 

Small drops of graphite black in an organic solvent (separation ~ 2mm) are carefully 
put on the sample and allowed to dry for about 30 minutes. 

2.1.4 Surface treatment 

arious treatments are done to change the surface and the changes in the surface stares 
are measured. These are light soaking, anneahng, native oxide etching, oxide growth 
and treatment with nitrogen and hydrogen plasma. 

Light soaking is done by white light (« 20 mW/cm^) from 300W tungsten-halogen 
lamp for 15 hours. A red filter is used to av'oid heating of the sample. Etching of the 
native oxide is done using HF (15% in deionised water) for 10 minutes. Oxide on the 
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surface is re grown by treating the sample with concentrated sulfuric acid, rinsing in 
deionised water and exposure to isopropyl alcohol, followed by annealing (ISO^C, 2h). 
Finally, we also treated the a-Si:H sample with nitrogen plasma (300K, 40mVF/cm2,dc 
plasma for 2 hours) and hydrogen plasma (300K, 40mW/cm^dc plasma for 2 hours). 

2.2 Characterisation of a-Si:H Films 

Films are characterised structurally by x-ray, optically IR and visible-UV spectroscopy 
and electrically by dark conductivity (cr) as a function of temperature and voltage, room 
temperature photoconductivity (cTp/i) and thermopower as a fimction of temperature. 

2.2.1 Structural Characterisation 

a-Si:H films are characterised for their amorphous nature using x-ray diffraction. The 
x-ray diffractogram gives a broad peak, characteristic of amorphous material. Electron 
diffraction consists of diffused rings, which also shows that the films are amorphous. 

2.2.2 Electrical Characterisation 

All doped and undoped films are characterised by dark and photoconductivity (cr and 
(Tp/i) at room temperature. For o-pf, measurement, light from a 300W tungsten halogen 
lamp (80V, 4amps) with a red filter is used. It gives about 30ml47 cm^ of intensity on 
the sample. We find a ^ and cTph ^ for our good 

quality undoped a-Si:H. 

Measurement of activation energy (E^) Activation energy of the dark conductivity 
is obtained by measuring cr as a function of T between 300 K and 450 K in vacuum. The 
details of the measurement chamber are described elsewhere [6]. Fig. 2.2 shows log cr 
versus 1/T curve for undoped a-Si:H thin film (# XI). is found to be « 0.75eV. 
Measurement of thermopower (S): For the mecisurement of thermopower (5). 
coplanar nichrome electrodes, about 1cm long and 0.5 cm apart, are evaporated on the 
top surface of the film. The sample is placed on two separate copper blocks, containing 
two independent heaters. The temperatures of the two ends (Ti and T 2 ) are monitored 
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10^/T (K'^) 


Figure 2.2; log cr plotted as a function of 1000/T for an undoped a-Si:H film. The 
activation energy (Ecr) is (0.75±0.05)eV. 
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using two copper-constatan thermocouples, placed on each electrode. By varying the 
heater currents, a temperature difference AT (= Ti - T 2 ) varying between 0 and ± 
5K could be achieved. Corresponding thermo emf (AV) is measured. Slope of the 
plot of AV as a function of AT gives S. These measurements are done in the partial 
pressure (0.5 mbar) of helium, for better thermal contact between the sample and the 
thermocouple. More details are given in Ref [6]. The thermopower is found to be 
negative and its magnitude increases with the temperature. These observations are in 
agreement with others [9]. Fig. 2.3 shows S as a function of 1/T for a P doped a-Si:H 
sample (#5). The slope is Ej = 0.13±0.03ey. 

2.2.3 Transmission Measurements 

Optical transmission (400-1200nm) measurements are used to obtain refractive index 
(n(A)), thickness(d), optical gap (Eg) and absorption coefficient (Qr(A)) of a-Si:H de- 
posited on 7059 glass. An accurate determination of a;(A), d and Eg using transmission 
mecLsurements is very important for our work because a(A) values are used for fixing 
the a scale in the CPM. Also, d is used in obtaining a from PDS signal and Eg values 
are used in the DOS modeling. In the present work we have used the transmission 
mecisurements and analysis followed by Swanepoel [10], which gives an accuracy of 1% 
for films with good thickness uniformity as claimed by the author [10]. The analysis 
and formula.e used, are given in .Appendix A in a closed form. 

Transmission spectrum is measured using a Hitachi - 2000 spectrophotometer for 
400nm < A < 900nm and a home made setup for 400nm < A < 1200nm. Similar values 
of transmission are obtained from the two set ups. Fig. 2.4 shows the home made set up 
used. It consists of a 75W quartz halogen lamp and a grating monochromator (Bausch 
& Lomb). The monochromatic beam is focused onto the sample using a lens. A high 
pass order sorting filter (Hoya-64) is used to remove the second order radiation from the 
monochromator output. To decrease the noise from stray light, we use a chopper (SR- 
540) and a lock-in amplifier (SR-510). The sample is mounted on a fixed, opaque plate 
with a hole of ~ 5mm diameter. The transmitted intensity is measured using calibrated 
photodiodes as shown in Fig. 2.4. We use a Si photodiode for 400nm < A < 1050nm 
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and a Ge photodiode for A > lOSOnm. After measuring the transmitted intensity' 
through the sample, the sample is removed from the plate carefully without disturbing 
the setup. Intensity of the incident light is now measured. 

The set up is fully automated. The automation consists of reading the lock-in amp- 
lifier data using a Personal Computer (PC) with serial port (RS 232C) communication. 
Also, a stepper motor controlled by the same PC is used to drive the monochromator 
grating as shown in Fig. 2.4. Details of the circuit used for controlling the stepper 
motor are given in Appendix B. The computer programme for the automation is given 
in Appendix C. Fig. 2.5 shows [a.huY^'^ as a function of hu for an undoped a-Si:H film 
(No. 5/94). The Tauc’s optical gap [11] is found to be (1.75 ± O.lOeV). The thickness 
calculated from the interference fringes is found to be 1080±50 nm, and is confirmed 
by a telestep apparatus. 

2.3 Sub- Gap Absorption Measurements 

The lowest a for 1 /^m film that can be measured using our transmission set up described 
above is about 50 cm~T However, a in the sub-gap region can be smaller by a factor 
of 100 or smaller. Hence more sensitive techniques are required. We have measured 
sub-gap a using two techniques namely photothermal deflection spectroscopy (PDS) 
[12,13]. and constant photocurrent measurement (CPM) [14,15]. 

2.3.1 Photothermal Deflection Spectroscopy (PDS) 

2.3. 1.1 Principle 

The PDS technique is based on the detection of heat produced when a sample absorbs 
sub band gap light called the pump light. As shown in Fig. 2.6, the pump beam 
(marked 1) is focused on to the sample, kept in carbon tetra chloride (CCI4). As the 
sample absorbs the pump light, the probe beam (marked 2 in Fig. 2.6) is deflected 
by the thermally induced changes in the refractive index of CCI4, which is in thermal 
contact with the sample. The deflection of the probe beam is mea.sured using a position 
sensitive detector (PSD). The deflection of the probe beam results in a signal (S) in the 
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> Light 

► Electrical signal 

Figure 2.4: Setup for Transmission measurements in visible-near IR region. Unfilled 
arrows indicate the light rays, whereas, the filled arrows indicate the electrical signal. 
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PSD and associated electronic circuit. The signal can be expressed as [12]: 


S = - exp{-ad)] (2.1) 

where, d is the film thickness and So is the PDS signal in the saturation region, since S = 
So when ad » 1. Thus, knowing So from the saturated region of the PDS signal and d 
from the interference fringes in the transmission measurements [10], a can be obtained 
using Eq. (2.1). However, this leads to a larger uncertainty in the determination of 
absolute value of a. This can be improved by matching the value of a in a wavelength 
region in which the transmission measurements can also give absolute value of o;, with 
greater precision. Details of the experimental set up [13] used in our laboratory are 
given below. 


2.3. 1.2 Instrumentation 

Fig. 2.6 shows the schematic diagram of the PDS setup fabricated by us. It consists of a 
monochromatic sub band gap light (pump beam 1 in Fig. 2.6), a He-Ne laser probe light 
(probe beam 2 in Fig. 2.6), a Position Sensitive Detector (PSD), an electronic circuit 
to enhance the signal to noise ratio and a lock-in amplifier to measure the signal. The 
pump light is from a high intensity grating monochromator which uses a 75W quartz 
halogen lamp as a source. It is chopped at a low frequency (9Hz in our case) to reduce 
the thermal lag. The monochromatic beam is focused on to the sample kept in a quartz 
cell containing pure and triple distilled CCI4. A long pass filter is used to filter out the 
second order from the grating monochromator. The probe beam (Fig. 2.6) is from a 
low power (2mW) He-Ne laser with good pointing stability (ours is better than 0.03 m 
rad/hour). The probe beam is focused using a short focal length (~ 4 cm.) lens near 
the center of the sample, and is allowed to pass nearly grazing the sample. The PSD is 
a quadrant photodiode. Since we need to detect the deflection of the probe beam only 
along one direction, two diodes out of the four in the PSD (which are perpendicu. ar 
to this direction) are shorted together to get two diode pairs. Details of the mounting 
of the PSD and the equivalent circuit diagram of the PSD are given in Appendix D. 
Actual components used by us are shown in Fig. 2.6. 
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The difference of the two outputs from the PSD is proportional to the PDS signal. 
These signals (called A and B) are processed in an analog electronic circuit to get (A- 
B)/(A+B) and amplified. Details of the divider circuit are in Appendix E. The signal 
is extremely small, and utmost care is required for good performance. The following 
points were kept in mind while designing our PDS system: 

The deflection medium (liquid) should be such that it does not react with the mater- 
ial under investigation. Also, it should not absorb light up to near IR and should have 
a large rate of change of refractive index with temperature (dn/dT). For a-Si:H films. 
ecu satisfies all these conditions and has dn/dT = 8 x 10“‘‘/K. However, special care 
has to be taken to handle CCI4, because it is highly volatile and can be a health hazard 
upon inhalation. As an alternate, Fluorinert FC72 (3M Product) can be used. This 
Hquid is less toxic, but has about a factor of 2 lower dn/dT (dn/dT = 4 x 10~'*/K) 
compared to CCI4 and also induces considerable change in the phases of lock-in signal 
[16] for photon energies < 1.4eV. Therefore, we have used CCI4 in a well ventilated 
room. Further, dust particles in the deflecting medium can cause tremendous noise. It 
is found that triple distillation of CCU a-nd clean environment while handling it helps. 

Finally, the mechanical vibrations can be a great source of noise. These are sup- 
pressed by keeping the experimental setup on a vibration free table having good atten- 
uation characteristics, particularly at low frequencies (10-20 Hz) (we used a vibration 
free table from Newport, Model M-MST-48-8). 

2.3. 1.3 Automation of PDS 

The PDS signal is very low and noisy. Hence long integration times (10-30 sec.) are 
required to reduce the noise. Also, to get a reasonable data set. one has to take readings 
at close intervals of wavelength. This is very tiresome and automation is desirable. 

Our automated PDS set up is shown in Fig. 2.6. The main features are essentiaOy 
the same as those m the manual mode described in Section 2.3.I.2. Additional arrange- 
ments for the automated measurements have been made using a computer (PC - 386). 
The PC has the following functions: 

a) Initialisation of the monochromator at the desired initial wavelength by moving the 
monochromator with the help of a stepper motor, run by the PC. 
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b) Reading tbe PDS signal from the lock-in amplifier (SR 510) using the serial port of 
the PC and RS 232 C communication. 

c) Moving the monochromator by desired steps using the stepper motor. 

The details of the circuit and software for running the stepper motor are given in 
Appendix B. 

The computer programme for the automation (Appendix C) consists of the following 
subroutines: 

i) Subroutine to read the data through Serial Port (for the lock-in amplifier SR 510). 

ii) Subroutines to run the motor driver circuits. 

2.3.2 Constant Photocurrent Measurement (CPM) 

2. 3. 2.1 Principle 

The real problem, when we try to estimate a from the photocurrent measurements, is 
that of changing life time of the electrons (r) with the changing Fermi level split [17]. 
In CPM, however, the photocurrent is kept constant by changing the intensity of light 
as the energy of incident light {hu) is varied in the sub band gap region. Therefore, 
the same set of recombination centers govern the free carrier lifetime and hence r is 
constant. The intensity that gives the predetermined constant value of the photo current, 
can be related to a as discussed below: 

For imdoped a-Si:H, the photocurrent is mainly contributed by the electron current 
because the fj.T product is about 10 times larger for electrons than holes [17]. The 
photocurrent density (J„) is given by [18] 

Jn = e np/i(l — R)[l — exp{—a d)] 77 pr Ff L. (2.2) 

Where e is the electronic charge, the number of incident photons per unit area 
per unit time, R is the reflection coefficient, d the film thickness, rj the photogeneration 
efficiency, ^ the electron mobility, r the mean lifetime of the electrons (recombination 
lifetime), F the applied electric field and L is the separation between the electrodes. 

In the sub gap region R is a slowly varying function and can be assumed to be a 
constant, pr product depends upon the position of quasi Fermi level, and is a constant 
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in this case. The value of rj is found to be « 1 at 300K [17]. However, by the comparison 
of CPM derived a with PDS derived a, it is found that r]{hu) does not decrease even 
in the sub gap region (with an uncertainty of about a factor of 2) [17]. Thus we obtain 
from Eq. (2.2) 


Jn = constant x [1 — exp{—ad)] Up/, 


(2.3) 


Now, for sub band gap region ord n 1, hence [l-exp{-ad)] reduces to ad. Thus for 
constant photocurrent Eq. (2.3) reduces to 


, , , constant 

' ( 2 -^) 
Therefore, by measuring the intensity that keeps the photocurrent constant at a 
predetermined value, we can find the spectral dependence of a up to a constant using 
Eq. (2.4). This constant is obtained by matching the CPM curve with the a obtained 
by the transmission in the region of higher energy where ad Ri 1. 

2.3. 2. 2 Instrumentation 

Fig. 2.7 shows the setup for CPM. Here the light source is the same as the pump beam 
part of the PDS apparatus as described already in the Section 2.3.I.2. A beam splitter 
(10 /o reflection) is used to monitor the intensity of the incoming beam using a calibrated 
photodiode (we used a Si photodiode, in the energy range 1.8eV > hu > l.UV and a 
Ge photodiode, in the energy range L2eV > hv > O.Sel/ ). The constant value of the 
photocurrent is read across a lAffi resistor using a lock-in amplifier. A bias of lOV 
( 'ectrode separation 0.2cm.) is used for the measurement of photocurrent (typically 
10 ^°.4). Actual components used by us are indicated in Fig. 2.7. 

2.3. 2. 3 Automation of CPM 

In CPM, photocurrent is kept constant by controUing the power to the monochromator 
lamp in a feedback arrangement at every wavelength of measurement. This reduces 
labour and gives better accuracy, compared to the manual measurements. 
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Figure 2.7: The setup for Constant Photocurrent Method. SMI and SM2 are the 
stepper motors. Lock-in 1 is an analog lock-in amplifier and hence, an A/D card (PCL 
812) is used to read this signal through PC. Filled arrows are the light rays, whereas, 
the unfilled arrows represent the electrical signal. 










Our automated CPM setup is shown in Fig. 2.7. The main measurement sections are 
essentially the same as in the manual mode described in Section 2. 3. 2. 2. The following 
are the additional arrangements for the automated measurements: 

a) Initialisation of the monochromator using the stepper motor ( SM2 in Fig. 2.7) used 
for driving the monochromator (Appendix B). 

b) Reading the photocurrent across IM Ohm resistance in the lock-in amplifier [lock-in 
(1) in Fig. 2.7] using A/D card and the PC. 

c) Comparing this photocurrent with the predetermined constant value of the photo- 
current, and moving the variac attached to the monochromator lamp using the stepper 
motor (SMI, Fig. 2.7) one step in the desired direction and comparing the new pho- 
tocurrent with the predetermined value. Repeating this process step by step till the 
photocurrent in the sample is within 2% of the predetermined photocurrent value. 

d) After reaching the predetermined photocurrent value, the intensity of the lamp is 
measured using calibrated photodiode and lock-in amplifier (lock-in 2 in Fig. 2.7). 

e) The monochromator is moved to another wavelength using the stepper motor (SM2). 

The movement of stepper motor is further discussed in Appendix B. The computer 

program for the automated CPM measurements (Appendix F) has the following sub- 
routines: 

>) Subroutine to read the data through Serial Port (for the lock-in amplifier SR 510). 
ii) Subroutine to read the data throgh A/D card (for the analogue Lock-in Amplifier). 

Subroutines to run the motor driver circuits. 

2.3.3 Dual Beam Constant Photocurrent Measurements (DB- 
CPM) 

CPM IS sensitive only to those optical transitions which end in the conduction band 
'Wended states. But for the determination of correlation energy (U), informs.tion about 
' optical transition ending in the localised states is also required. This .s achieved 
- changing the occupancy of the localised states using a dc bias light. This can help 
'« obtaining U and the shape of the deep gap states, as discussed later in Section 3.3. 
Dual beam constant photocurrent measurements (DBCPM) [19] is ac CPM done 
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in the presence of a highly absorbing dc light, known as the pump light. We have 
used another high intensity monochromator as the source for the intense dc light. The 
details of the instrumentation and automation of the set up are the same as discussed in 
section 2.3.2. We have used a quartz halogen lamp (120 W) along with a high intensity 
grating monochromator (Oriel), as the pump light. The wavelength used is 650nm. 
The intensity of the pump light is approximately 20 times that of the chopped (probe) 
light. Fig. 2.7 also represents the DBCPM set up without the bias light. 
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Chapter 3 


Bulk States 


Sub-gap absorption [1-4] is a fundamental process that can measure the density and 
the energy position of all the states in the mobility gap, provided the matrix elements 
of the optical transitions are known. These are usually assumed to be constant and 
independent of photon energy (hi/). In a-Si:H, disorder gives rise to the Urbach edge 
[5], whereas, other defects like the dangling bonds give deep states, responsible for 
the absorption at smaller hw. The a associated with the absorption in the deep gap 
states is too small to be measured by conventional transmission method. For this, two 
techniques namely, constant photocurrent measurement (CPM) [1,2] and photothermal 
deflection spectroscopy (PDS) [3,4] are employed. Details of CPM and PDS set ups 
along with the procedure to obtain ct{hi/'^ from CPM and PDS spectra are given in 
Chapter 2. 

Sub-gap absorption coefficient measurements have been used by earlier work- 

ers to obtain the density of bulk states in a part of the mobility gap of a-Si:H, using 
CPM [1,2] and PDS [3,4]. However, states in the entire mobility gap using sub-gap 
absorption have not been obtained. In particular, the slope of conduction band Urbach 
edge has not been measured. Further, there is a large variation in the value of correl- 
ation energy (U)(0<U<0.5eV) [6-8] measured using various techniques, although, it is 
generally agreed that U is not negative [6]. 

Estimation of DOS using a(hi/) requires the knowledge of the shape of DOS a priori. 
But different authors [9,10] prefer different shapes for the mid-gap states. In particular. 
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Gaussian [9] and exponential [10] shapes for the mid-gap states have been proposed. It 
is noted that using the sub-gap absorption measurements one can not obtain the shape 
of the mid-gap states uniquely [11,12]. 

We have used a^pm and Oipds data to obtain the density and distribution of the bulk 
states in the entire mobility gap for lithium doped a-Si:H [a-Si:H(Li)] and phosphorous 
doped a-Si:H [a-Si:H(P)]. Lithium (Li) is an interstitial n type dopant in a-Si:H [13], 
unlike phosphorus (P) which is a substitutional. Thus, it is interesting to see whether 
dopings by P and by Li change the DOS in a similar manner. Also, it is worthwhile 
to see if the uncertainties that plague the DOS in a-Si:H(P) remain for a-Si:H(Li) as 
well. One such uncertainty is the position of D~ in a-Si:H(P). Street et al [9], based 
on luminescence and transport measurements, conclude from their ’’Basic Model of 
DOS” that the positions of the D" peak below the conduction band edge, in a-Si:H 
and a-Si:H(P) are the same. On the other hand, Kocka et al [14], based on sub-gap 
absorption measurements, and are of the view that on doping, D~ moves towards the 
center of the mobility gap, because of intimate bond pair formation between dangling 
bonds and dopants. This is known as “Intimate Pair Bond Model” [14]. However. Li 
may or may not form an intimate pair with the dangling bond, as P does. Therefore the 
position of D~ peak in a-Si:H(Li) may be different from that in a-Si:H(P). A comparison 
of the two different dopings shows that these are similar as far as the density of bulk 
states are concerned. Further, we find that even after a careful deconvolution of sub- 
gap absorption data and detailed error analysis in the case of a-Si:H(Li), we can not 
distinguish between the models by Street et al [9] and Kocka [14], regarding the position 
of D . We try to remove this uncertainty by combining the sub-gap absorption data 
with the conductivity activation energy data. 

In Section 3.1, we describe the principle and analysis of sub-gap absorption data 
to obtain DOS in a-Si:H. First we describe various shapes of the DOS distribution 
proposed in the literature and sub-gap transitions. Determination of slopes of Urbach 
edjjOs and deep gap states are discussed then. In the next subsection, we give the 
analysis of DBCPM results for the determination of U and the shape of deep gap states. 
In Section 3.2, we study the effect of doping on the density and distribution of bulk 
Results on P and Li doped a-Si:H with various concentrations of dopants are 
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given. Section 3.3 gives the summary and conclusions of this Chapter. 

3.1 Undoped a-Si:H 

3.1.1 Principle and Analysis 

3. 1.1.1 DOS distribution model and sub-gap transitions 

Although, it is generally agreed that the states near the band edges are exponential, 
various authors prefer dilferent shapes for the deep gap states [9,10]. Two shapes namely. 
Gaussian [9] and exponential [10], have been proposed. The schematics of these DOS 
distributions are shown in Figs. 3.1 and 3.2. We discuss them one by one. 

Fig. 3.1 shows a model of DOS distribution [9] for undoped a-Si:H proposed to ex- 
plain many measurements including the sub-gap absorption in the literature. Also shown 
are the various possible optical transitions (marked T\ to Tz). The model DOS has an 
exponential shallow region (Urbach energy region) and deep defect region (dangling 
bond states). The neutral dangling bonds {D°) are separated from the doubly occupied 
dangling bond states {D~) by correlation energy AE. D° and D~ are both assumed to 
have the same shape and number. The Fermi level is at the minimum of the DOS, as 
shown in Fig. 3.1. Fig. 3.2 shows another model of DOS, with the difference from the 
first that the deep gap states in this case, are exponential with a slope different from 
those in the Urbach edge region. Similar optical transitions (Ti to T5) are shown in 
this case also. 

In undoped a-Si:H, the photocurrent at room temperature is dominated by electrons. 
Therefore, CPM will p;'imarily see only those optical transitions which end in the con- 
duction band extended states, nameljs transitions Ti, T2 and T5. The transitions T3 
and T'4, on the other hand, take the electrons to localised final states and hence are not 
seen by CPM. All these transitions, however, generate heat and will clearly contribute 
to a measured by PDS. Here, transitions from one localised state to another localised 
state are neglected because the matrix elements for these transitions are likely to be 
small. Further, CPM is insensitive to the surface states because the electron hole pairs 
generated at the surface quickly recombine without contributing to the photoconduct- 
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Figure 3.1; Distribution of density of bulk states[after reference 9] used in the model cal- 
culations along with the optical transitions in a-Si:H. It consists of exponential Urbach 
edges and Gaussian dangling bond states. PDS measurements see all the transitions Tj 
to Ts whereas CPM ignores the transitions T 3 and T 4 . E/ and Ei have been measured 
from Ec- 






ivity [15,16]. Hence, for the transitions in Figs. 3.1 and 3.2, we can write in an obvious 
notation 


apj, = a(Ti} + a(T 2 ) + 0(^3) + 0(74) + a(r5) + a(Ts) 

and, 

0!cpm = ^C^l) + <^(^2) + CulTs), 

here, r, refers to the transitions involving surface states (not shown in Fig. 3,1 and 

■3.2). Absorption attributed to the transitions T, will be used in Chapter 4 for the 

determination of density of surface states in a-SiiH. Thus, various regions of absorption 

in CPM and PDS spectra shown in Fig. 3.3 may be represented by the following 
transitions: 

Region I: 


Region II: 

^pds — ^cpm — Qf(Ti). 

(3.3) 

and. 

a^^, = a(T 2 )^a(T 3 ), 

(3.4) 

Region III 

^cpm ~ 01(1^2) • 

(3.5) 

and, 

(^pds = 01(74) + 01(75) + a(Ts), 

(3.6) 


^cpm = Q'(J 5 ). 

(3.7) 

3.1.1.2 Determination of Slopes of Urbach edges 

In region 11 (Fig. 3.3). Itnown as the Urbach edge region. .s attributed to the 

nd on band eictended states (T, in Figs. 3.1 and 3.2). The absorption curve can 
oe titted to an exponential curve [5] 


^urbach — O’^ expikul 

(3.S) 
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Figure 3.3: log a as a function of iii^ as measured by CPM (filled circles), PDS (open 
circles) and Transmission (crosses) for an undoped a-Si:H sample. Peaks in the curves 
arise from interference effects. 
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where, Eot, is the slope of valance band Urbach edge. In this region, is more than 
Ocpm- This is because of the additional transitions from VB extended states to the 
conduction band Urbach edge localised states (T 3 , Figs. 3.1 and 3.2). It is clear from 
Eqs. (3.4) and (3.5) that, if the plot of log{apds - ^cpm) as a function of hu is a straight 
line (corresponding to the transition T3), its slope gives Eoc- 

3. 1.1.3 Determination of Deep Gap States 

a. Gaussian shape of mid-gap states (Fig. 3.1) 

In the one electron approximation, the absorption coefficient a(hu) for transitions 
from initial states gi to the final states p/ is given by the relation [17]: 

a(H = ^ 9fi^ + - /(e + hu)]de. (3.9) 

The integration is carried over all the initial states p,(e) and final states gf(e + hu) 
separated by the energy hu. f is the Fermi Dirac function. We take zero temperature 
statistics and hence, the states above the Fermi level are taken to be empty (/(c) = 
0 ) and the states below the Fermi level are filled (/(e) = 1 ). The factor C contains 
the spectral dependence of dipole matrix element which have been shown [18] to be a 
weak (within a factor of 2 ) function of hu for transitions from localised states as well as 
extended states. C is taken to be a constant for transitions between the localised states 
to the extended states and also for the transitions between the extended states to the 
extended states. The constant C can be estimated by considering the transitions from 
the extended valance band states to extended conduction band states and assuming 
these states to be parabolic. It is further assumed that g(e) at the band edges is given 
by free electron density of states (6.7 x lO^^cm-^eU-^) [ 2 ]. Hence. 

P.-(e) = 6.7 X 10^' (e)('/ 2 ) = 

where, e is the difference between the energy of the states from tHe band edge. So, for 
the band to band transitions Eq. ( 3 . 9 ) gives [ 12 ] 

a{hu) = —(6.7 X 10^^)^^ ■ eU2 _ ^^ 1/2 
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Thus measurement of a for hu > Eopt gives an estimate of C using Eq. (3.11). It 
is evident from the Eq. (3.11) that q is a convolution of all the initial states and the 
final states, which obey the conservation of energy. So, to get information about DOS 
in the mobility gap we have to assume the shape of DOS a priori. The plateau region 
in a(hu) curve (region III in Fig. 3.3) is because of the transitions from Gaussian D° 
level to the CB extended states (transition Tg) in Fig. 3.1. Under this assumption, p,- 
is given by the relation: 


gi{e) = N, exp[-~{Ei ~ hiyf/2W% (3.12) 

where is the height of the Gaussian peak, Ej is the position of the peak below the 
conduction band edge and W is the half width of the distribution. Using this distribution 
of DOS, Eq. (3.9) reduces to: 

Oinff{hu) = ^ 6.7 X lO^UV^ + ]e^/^de, (3.13) 

where integration limits (a,b) is taken from 0 to a reasonable value 2eV [1] and adifj{hu) 
= a{hi') - extrapolated value of o: from region II (Fig. 3.3). is calculated using 

Eq. (3.13) in the region III (Fig. 3.3), taking g{Ej), W and Ei as parameters. The 
best fit between the experimental and calculated values gives the fit parameters, 
b. Exponential shape of mid-gap states (Fig. 3.2) 

For the exponential distribution of the mid-gap states (Fig. 3.2), g,(e) is given as 

gi{e) = ,Vea;p(-^^g— ), (3.14) 

where, N is the extrapolated value of the mid-gap states at the valance band mobility' 
edge and E^o is the slope of deep states. For this distribution of the deep states Eq. 
(3.11) becomes 


aaiiihi,) = ^(6.7 x 10^‘)iV (3.15) 

The parameters N, Eod and Ej^ obtained from the best fit between the exper- 
imental (Xdiff and the calculated a*//. 
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3.1. 1.4 Dual Beam CPM Measurements 


a. Principle and theory 

Dual beam CPM [6] is ac CPM (using chopped light) in the presence of a strong 
dc (unchopped) light. The weak light from a monochromator, chopped at a frequency 
between 5Hz and 333Hz, is the probe beam. The strong dc light is the pump beam. The 
details of the experimental set up are given in Section 2.5. Fig. 3.4 shows normalised 
spectra measured by dc CPM (curve 1), and DBCPM with the weak light chopped 
at 333Hz (curve 2), 33Hz (curve 3), 9Hz (curve 4) and 5Hz (curve 5) respectivelv for 
photon energies O.SeV to l.SeV. All these spectra are normalised to the same value of 
a in the region of photon energy 1.8 < hu < 1.6el/, obtained from the transmission 
measurements [19]. 

DBCPM spectra (Fig. 3.4) has two interesting features. First, the spectra at oHz 
chopping frequency has a prominent hump at low photon energy (hi/ < l.OeV), which 
decreases with increasing chopping frequency. Second, there is a strong dependence of 
the entire DBCPM spectrum on the chopping frequency. 

We first explain the hump in DBCPM spectra at low chopping frequency. I'pon 
measuring the dual beam spectrum using PDS, on the same sample, at the same fre- 
quency, we find that the hump is absent (DBPDS data not shown) [6], This indicates 
that the hump may not be true absorption and may arise because of a change in ac 
photoconductivity in the presence of dc light [6]. Therefore, we first discuss the dual 
beam photoconductivity in a-Si:H. It is assumed that the pump beam is so intense that 
the defect states lie within the quasi Fermi levels, and Ej, as shown in Fig. 3.-5. 

In DBCPM, the excess electrons in the conduction band generated by the chopped 
light, in the presence of the ])ump beam, are given by [6] 


n-GT = aFr, 

where, r is the recombinatroa life time for the strong d.c. pump beam alone and G is 

the generation rate from the pump light and F is the intensitt’ of the weak chopped 
light. Differentiating Eq. (3.16), we get: 
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hi^(«V) 


Figure 3.4; Absorption coefficient a as a function of )ii/ obtained by dc CPM (curve 1), 
DBCPM with the weak (probe) light chopped at 333Hz (curve 2), 33Hz (curve 3), 9Hz 
(curve 4) and 5Hz (curve 5), for an undoped a-Si:H (Sample No. 18). 
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log [n(E)] 



Figure 3.5: A schematic of density of localised states in undoped a-Si:H along with the 
optical transitions in the dual beam configuration. E /„ and E /p are electron and hole 
quasi Fermi levels respectively (after reference 6). 
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An = GAt + aAFr, (3-17) 

where, At is the change in the recombination life time due to the weak chopped light, 
AF is the intensity of the weak chopped light. The hump in DBCPM spectra at low 
frequencies of chopping can be imderstood if At has two components. We write 

At = Ati + At2 . (3.18) 

We discuss Atx first. At^ is identified with the known power law of photocurrent 
{Iph.(x F'^), with 7 < 1, for a-Si;H [20]. Therefore for An part, 

(3.19) 

(3.20) 

An = GAt 2 + jarAF. (3.21) 

If At 2 = 0, Eq. (3.21) correctly yields a, as in CPM. But the observation of hump 
suggests that At 2 7^ 0. We postulate that At 2 results from a change in the occupation of 
recombination centers through excitations by the probe light between the defect states 
and the conduction and valance band. The hump, therefore, requires two additional 
transitions in DBCPM, which are insignificant in single beam CPM. These transitions 
are identified as follows: 

Fig. 3.5 shows schematic density of states along with the optical transitions (Ti to 
T5), in the dual beam configuration. For this discussion we have chosen the Gaussian 
distribution of deep states. Transitions T4 (Fig. 3.5), with a lower onset energy convert 
some of the D"*" states to D°. This increases the photocurrent because, for electrons. 
D+ is about 10 times more efficient recombination centers [18] than D°. Transitions 
T2 with a higher onset energy, on the other hand, change some of the D° to D"*" and 
decrease the photocurrent. This explains the hump in the DBCPM spectrum. 
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b. Determination of correlation energy (U) 

The correlation energy (U) can be obtained by analysing the DBCPM spectra at 
high chopping frequency (333 Hz). At this frequency of chopping, the second term 
in Eq. (3.17) dominates and the first term is insignificant, because the hump in the 
DBCPM spectrum disappears. DBCPM spectrum at 333Hz arises from the transitions 
from the filled D“ states to the conduction band extended states (T5, Fig 35 ) 
fit between the experimental and calculated a{hu) gives the characteristic energy (E^) 
of the transitions T 5 . This gives U = F, - where E, is the characteristic energy 
for the transitions T 2 (Fig. 3.5) and has been obtained from the analysis of dc CPM 
spectrum. 

c. Shape of mid-gap States 

In the above calculations we have used Gaussian distribution of DOS [ 9 ] (Fig 3 5 ) 
However, another shape of the DOS (exponential) (Fig. 3 . 2 ) has also been suggested 
[10], We use the two DOS to fit dc CPM and DBCPM data to see which shape fits 
our data; The parameters Wo, W and B. (Fig. 3.5) are varied to git, best fit between 
the experimental ar.-„(A,.)(= p i„ region 11 - a,exp(-/u./f;,) extrapolated in region 11 
(Fig. 3.4)1 the calculated 0 * 7 /, for dc CPM curve. We get a good fit between the 
wo. Also, good fit is obtained for exponential shape of DOS. We are, therefore, unable 
to determine the distribution of DOS uniquely from the CPM measurements. A similar 
conclusion was reached by Payson and Cuba [ 11 ], for their PDS results. Although, a 


ussian distribution for the mid-gap states explains the luminescence, time of flight 

and junction measurements [9], it remains to be seen whether some other distribution 

of DOS also fits the observed results. 

We now fit DBCPM data at 5 Hz chopping frequency using the Gaussian and expo- 
nential shapes of DOS. The hump at low photon energy (Fig. 3.4) is attributed ,0 the 
transitions T, and T, (Fig. 3.5), as discussed in the Section 3.1.1.4a. The first term 
m Eq. (3.21) IS the dominant term for DBCPM spectra obtained for low frequency of 

chopping (5 Hz). An expression for the Aa (= a . at SHz ^ f ■ • tt p- 

^ ^dbcpm at oaz - ^cpm): m region II iFig. 

3.4) can be written as : 


Aa = -41 


%-nv -F(-eV2W=)* - exp(-eV2IV^)*], 


(3.22) 
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UJD) 



hi'(eV) 

Figure 3.6: The peaks in the experimental data (o o o) are from interference effects. 

These are averaged to obtain the mean experimental curve ( )- Curve marked 1 and 

2 are theoretically calculated using Gaussian and exponential deep states, respectivelv. 
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where, E. and E, are the characteristic energies for the transitions and respect- 
ively (Fig. 3.5). A is the constant of proportionality. The first term m Eq. (3.22) 
is proportional to the net number of D+ changing to D° and therefore, represents the 
increase in photocurrent because of the presence of pump light. Similarly, the second 
term is proportional to changing to D+, and represents the decrease in photocurrent. 
We already know Nq and W from the deconvolution of DC CPM spectrum. Hence, 
from the best fit between the theoretical and experimental Aa, we get Eb which is the 
characteristic energy of the transitions T 4 . It may be mentioned that Ea + Eb = Eg, 
the band gap of the material. Eg has been obtained from transmission measurements 
(Section 2.4) [19], also. This gives a check to the analysis. 

Aa can be calculated for the exponential shape of the deep gap states (Fig. 3.2), 
using the following relation : 

Aa = C{ exp{-e/EoD)dt- exp{-tl EoD)de], (3.23) 

where, C is the constant of proportionality. The Fermi level split is taken to be the same 
as in the case of Gaussian shape of deep gap states. This gives Ef — Ejp ^ 0.25eV. 
Hence, in this case there is no other fit parameter. 


3.1.2 Results and Discussion 


Slopes of Urbach edges 

Fig. 3.7 shows a[hu) measured by transmission (hz/ > 1 . 6 eV), CPM and PDS (in 
the region 11 and III). This figure is the same as Fig. 3.2, but the interference fringes 
have been averaged out. For all the CPM and PDS spectra from now onwards, we will 
show the fringes avei-aged out. For LSeK > hu > lAoeV (region II) the CPM curve 
can be fitted to an Urbach tail with Eov = 52±5meV [Eq. (3.10)]. Fig. 3.8 shows 
^^ai^pdsihu) - acpm.{^.u)) plotted as a function of hu in region II. Eoc is found to be 28 
AomeV. 


Mid-gap states 

adiff{hu) = oicpm{hu) — aoexp[—hu j Eov) in the region III of Fig. 3.7, is calculated 
using the method discussed in Section 3. 1 . 1. 3 for various choices of the Gaussian deep 


KAHrUt 




Tiliil 
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ho> (eV) 

Figure 3.7; a{hu) measured using CPM and PDS, for the same sample shown m Fig. 
3.3. The interference fringes in the cur\-es have been averaged out. 
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gap states parameters No, W and E,. The best fit between the experimental and 
calculated ns.;/ gives values of No, W and C, calculated using Eq. (3.11), is fomid 
to be 2x 10“^**. This compares favorably with the reported value [1,2]. Here is 
taken to be l.SeV as calculated by the Tauc’s plot (21) (Section 2,3.3). The following 
parameters are obtained for Gaussian deep gap states (Fig. .i.l), 

No = (1.8 ± 0.5) X 

El = (1.0±0.05)eK 


W = (0.14 ± 0.01)eV 

For the exponeatial shape of deep gap states (Fig.3.2), a satisfactory fit between 
calculated and experimental adi/f is found for the following DOS paiameters. 

iV = 4.0 X 10*®cm-'^eV'-^ 

EoD = 0.3eV, 

and, 

Epn — 0.6eV. 

adiff calculated using the two shapes of mid-gap states along with the experimental 
points are plotted in Fig. 3.8. 

Correlation Energy (U) 

Section 3. 1.1.4 gives how analysis of DC CPM and DBCPM wfith ac light chopped 
at 333Hz can be used to get characteristic energies Ea and E,i for transitions T 2 and T ,5 
respectively (see Fig. 3.5). We find E^ and Ei to be 0.10 ± 0.03eV and 0.85 ± 0.03eV 
respectively. This gives correlation energy U = 1.0 - 0.85 = 0.15 ± O.OfieV. 
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Shape of deep gap states 

As discussed above both Gaussian and exponential mid-gap states fit a(An) data 
satisfactorily. We are, therefore, unable to determine the distribution of DOS from CPM 
measurements alone. As dtscussed in Section 3.1.1.4, An is calculated for Gaussian as 
rvell as exponential deep gap states using E,s. (3.18) and (3.19) respectively. Fig. 
.3 9 shows An as a function of hn obtained exper.mentally as well as calculated using 
Gaussian and exponential shapes of the deep defect states. For the Gaussian shape the 
fi, parameters are E. = (1.0±.03)eV and E. = (0.80±.03). This gives the band gap. 
£ _£ -|-Ej = l8± 0.06eV. which compares well with the values calculated for the 

sample using transmission measurement [19]. This gives confidence in our analysis. 
The absence of a good fit, between the experimental and the theoretical values of Ac 
for the exponential shape of deep defect states leads us to believe that the hump in the 
DBCPM spectra can not be explained using the exponential shape of the deep delect 
states. Therefore, using a combination of DC CPM and the DBCPM spectra, we have 
been able to distinguish between the two shapes of the deep defect states and our data 
prefer the Gaussian shape. 


3.2 Effect of Doping on Bulk States 

3.2.1 Principle and Analysis 

3. 2. 1.1 DOS distribution and sub-gap Transitions 

Fig. 3.10 shows a schematic of the distribution of DOS for a-SiiH doped u type doped ^9] 
along with the optical transitions (Ti to T 4 ). Not shown in the figure are the transitions 
to and from the surface states (T^). Therefore, we can write as before, in the case of 
undoped a-Si;H (Section 3. 1.1.1): 

OiPDS = Q:(Ti) + oc{T2) + ociTz) + oc{T4) + oc{Ts) (3.-4) 


UCPM = a(r,) + a{T 2 ) + a(T,). (3.i5) 
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ha> 


Figure 3.9: Aa (= ccdbcpm at 5 Hz - acp. 
ergy. Solid circles are the experirnental points 
Gaussian and dotted line is the theoretical fit 







Figure 3.10: Distribution of density of bulk states used in the model calculations along 
with the optical transitions in n type doped a-Si:H. It consists of exponential valence 
band and conduction band Urbach edges with slopes and Eoc respectively. PDS 
measurements see all the transitions Ti to T 4 whereas CPM ignores the transitions Ts 
[after reference 9]. 
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3. 2. 1.2 Determination of Slopes of Urbach edges 

Slopes of valance band and conduction band Urbach, edges in doped a-Si:H are determ- 
ined using the procedure given in Section 3. 1.1. 2. 

3. 2. 1.3 Determination of Deep Gap States 

Deep gap states are determined by the deconvolution of aapmihv) for photon energies 
< l.SeV. The procedure is the same as discussed in Section 3. 1.1. 3 for undoped a-Si;H. 

3. 2. 1.4 Position of D” in a-Si:H(Li) 

In the "Basic Model” of Street et al [9], D~ in a-Si:H(P) is separated from D° in 
undoped a-Si;H by a correlation energy (U) of fsi 0.15eV. On the other hand, in the 
"Intimate Pair Bond Model” by Kocka et al [14] D° in undoped a-Si:H and D~ in a- 
Si:H(P) are located at the same position, because doping creates predominantly intimate 
pairs of Si-dangling bonds with dopants. The energy of these pairs lie deeper (measured 
from CB) than the energy of isolated Si dangling bonds. 

As discussed in Section 3.1.1, we can obtain the position of D“ states (Ej), W and 
N(E/) by deconvolution of acpmihu). .A.n error analysis is done for determining of E; 
accurately, as described below. Error (e) = (calculated value of a*// - measured 
value of a^iff at several points)'^ is plotted as a function of E/ for various values of W. 
N(E/) is kept constant. An absolute minimum in these curves gives the value of E/. 
In spite of careful analysis, we are not able to distinguish between the models of Street 
et al, and Kocka et al, because both models are within error bars. In order to resolve 
this issue, we measure the position of the Fermi level in these samples by measuring 
Eo-, the activation energy of dark conductivity. E^, can be taken to be approximately E. 
- Ef, within lOmeV or so [22]. The equality is not strictly valid because there may be 
a statistical shift of Fermi level cls well as other effects. It is clear from Fig. 3.10 that 
the DOS at the Fermi level g(Ef’) for a-Si;H(Li) can be written as 

giEp) = 2N,exp[-{Ei - E,)V{2W% (3.26 - 

It is seen from Eq. (3.22) that the %'alue of g(EF) strongly depends upon E/. We 
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Table 3.1: DOS parameters [Eov, Nj, Ej, W and Eoc) for a-Si:H(Li) with two Li 
concentrations and a-Si:H(P) with two P concentrations. 


Sample 

details 

E. 

(eV) 

(meV) 

Ni 

xlO"^® 

{an~^eV~^) 

El 

(eV) 

w 

(eV) 

Eoc 

(me Vi 

a-Si:H(Li)(#l) 

0.45 

70±5 

1.5±.l 

0.92±.06 

0.17±.005 

32±5 

a-Si:H(Li)(#2) 

0.38 

90±5 

3.1±.l 

0.95±.06 

0.18±.005 

34±.5 

a-Si;H(P)(#5) 

0.30 

85±5 

3.0±.l 

0.88±.06 

0.185±.005 

32±5 

a-Si:H(P)(#15) 

0.28 

90±5 

3.2±.l 

0.85±.06 

0-185±.005 

30±5 


calculate the value of g(EF) taking the value of E/ assuming the Street’s model and 
Kocka’s model to be correct. The value of g(EF) under preferred model will be in 
agreement with that obtained from other methods of determining gfEp-) [23]. 

3.2.2 Results and Discussion 

3. 2. 2.1 Lithium doped a-Si:H(Li) 

Figs. 3.11 and 3.12 show a{hp) for a-Si:H(Li) measured by CPM and PDS for samples 
#1 and #2 respectively. The interference fringes have been averaged out [22]. The 
sample #2 (E^ = 0.38eV) has a higher Li concentration than the sample #1 (E, = 
0.45eV), as is clear from the E^- value. 

Slopes of Urbach edges 

The slope of the valance band Urbach edge (E^v) is determined by Eq. (3.10) in 
region I (Figs. 3.11 and 3.12). E^^ is found to be (70 ± 5)me\’ and (90 ± 5)me\' for 
samples #1 and #2 respectively. 

g 3 13 shows the log ^ as a function of h.u in the region I for both the samples. 
Eoc IS obtained from the slope of this curve (see Section 3.2. 1.2). and found to be (32 ± 

o)meV and (34 ± 5)meV for the samples #1 and #2 respectively. The values of slopes 
of the Urbach edges are given in Table 3.1. 

Deep gap states 
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h x) (e V ) 

Figure 3.11: Absorption coefficier.t a a5 a function of hz/ as naeasured by CPM and PDS 
for a-Si:H(Li) (sample #1). CPM curve have been normalised by suitably matching 
CPM data with a{hi/) obtained from the transmission measurements in the region 1.8 
to 1.6 eV. The error bars shown in this figure are the same for Figs. 3.12 also. 
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Figure 3.13: adiff plotted as a function o 
angles). Theoretical fit ( — ) is based oi 
Please note the different scales for the two 






As discussed in section 3. 2.1. 3, acpm in the region II is anaivsed under thp aoc 

" '“Sumption 

of the Gaussian distribution of deep gap states. The Gaussian parameters g(^,) pj, 

and E/ are obtained from the best fit between the measured Q{hi/) and the theor t‘ 1 

a{hu). The values of these parameters for both the samples are given in Table 3 1 

As discussed in Seccion 2.1.2, we have made the Li dop,„g by vacuum evapoltioa of 

Li on uudoped a-Si:H. Although, exteuded annealing have been done to ensure uniform 

diffusion of Li throughout the bulk, secondary ion mass spectroscopy (SIMS) have not 

been done. So, the non uniform distribution of Li in the bulk can not be ruled out It 

is known that the the density of bulk states strongly depends upon the concentration 

of the dopant. Hence, the density of the bulk states obtained, shows an average value 
using the present analysis. 

3. 2. 2.2 Phosphorous doped a-Si:H 

Figs. 3.14 and 3.15 show o(An) for a-SiiH(P). sample #.5 and # 1-5 respectively meas- 
ured by CPM and PDS. Following the analysis discussed in section . 3 . 2.1 o slopes of 

> (Fi. 3.13 . 

f, I (Fig. 3 ., 5 )) plotted as a function 

of hn for samples #.5 and #15. The slopes give for the two samples. We «nd 

I E :rr V ^ 

_ (90± 5)meV and E„, = (30± 5)me\'. 

Qfcpm in region II (Fig 3 i-n 

, ^ parameters for the deep gap states, 

both the 

samples (#o and #15) along with the slopes of the Urbach edges are given in 


3.2.3 Position of D" in a-Si:H(Li) 


Fig. 3.17 shows o as a function of hn between 0 9eV and 1 3eV f Q. ir u 
Si:H(Li) with various concentrations of Li using CPM f " " 

tond parameters are obtained und ,h ’ ' ' ‘ 

DOS parameters are shown in Tal 

.J- The position of £>» in a-Si:H and £)- in 
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Figure 3.14: Absorption coefficient a as a function of hi/ as measured by CPM 
PDS for a-Si:H(P) (sample #5). 
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(e V ) 


Figure 3.15: Absorption coefficient 
PCS for a-Si:H(P) (sample #151 


function of hu 


as measured by CPM and 



a-Si : H(P) 



hVCeV) 

Figure 3.16: o;*-// plotted as a function of h.i/ for samples ^5 (circles) and ^15 (tri- 
angles). Theoretical fit ( — ) is based on Eq. 3. . Please note the different scales for 
the two samples. 



a-Si:H(Li) are obtained by a detailed error analysis. We find that prediction of the 
position of D~ in Li doped a-Si:H using the basic model of Street et al and intimate 
pair bond model of Kocka et al are within the error bars. However, g( Ep) for a-Si:H(Li) 
calculated using Eq. (3.22), is 2.5 xlO^' c'm~^eV~^ for the model by Street et al [9] and 
9.1 xlO^® cm~^eV~'^ for the model by Kocka et al [14]. We find that ^{Ef) in the case 
of Kocka’s model is unrealistically low, when compared with the values obtained by 
other methods for a similar (decided by E, value) a-Si;H(P) sample. Thus we conclude 
that the ’’Basic model” of Street et al [9] is probably tangible and is favoured by our 
data, for lithium doped a-Si:H. This looks reasonable because Li is an interstitial dopant 
m a-Si;H and hence may not form an intimate pair with the dangling bond states. 


3.3 Summary and Conclusions 

1. Sub band gap a has been measured using CPM and PDS on undoped and P and Li 

doped a-Si:H thin films between photon energies (hi/) 0.9eV and 2.5eV. is found 

to be larger than in the sub-gap region (hu < 1.8eV). This has been reported by 

others, as weU (24). The difference is explained by arguing that ,s not sens.tice 

to a part of the bulk states transitions and also the transitions to and from the surface 

states. Using CP.M and PDS for determining a as complementary techntques, we obtain 

density of bulk states in the entire band gap. The slope of conduction band Urbach 

edge has been obtained, which is not possible using CPM or PDS alone We have 

calculated the bulk states using a Gaussian as well as an exponeut.al mid-gap states 

Both gtve good fits to n(hv) data. The parameters obtained are in agreement with the 

terature. This shows that the shape can not be deternnned umquely using CPM and 
"Uo measurements. 

2. Shape of the mid-gap states, however, have been uniquely obtained us.ng DBCPM 

tneasurements at low chopping frequency Ha). A prominent hump appearrng in 

e spectrum be fitted using the Gaussian shape of the mid-gap states alone. So, 

we o .ude t at our data favour the Gauss, an shape of the mid-gap states. Further, 

DBCPM spectrum at 333Hx .s used to obtain the correlation energy (U) between the 

uvo electrons occupying the same defect site. We find U = 0 15 ± 0 06eV ' 

u.io ± u.UbeV, m agreement 
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Table 3.2: DOS parameters {Eov,Ni,Ei, and W), activation energy (E<,) and DOS at 
the Fermi level calculated using Eq. (4.5) assuming Street s model [^ 5 (£'f’)] and Kocka's 
model [gKiEp)] of DOS for a-Si:H (#18) and a-Si:H(Li) with two Li concentrations. 


Sample 

details 

Eoy 

(meV) 

W/xlO-^® 

Et 

(eV) 

w 

(eV) 

E, 

(eV) 

gsiEp) xlO-^^ 

gK{EF)xl0~^^ 

a-Si:H(#18) 

45 

5.85 

1.03 

0.119 

0.81 

- 

- 

a-Si:H(Li)(#l) 

70 

150 

0.92 

0.170 

0.45 

24 

0.90 

a-Si:H(Li)(#2) 

90 

310 

0.90 

0.180 

0.38 

28 

0.92 


with the Literature [6]. 

3. Bulk states in the entire mobility gap have been calculated using occPM and apQs 
measurements for P doped a-Si:H with two P concentrations. We find an increase in 
Eo„ with the increasing P concentration. This is also found by others [25]. Also, the 
density of deep gap states increases with increasing P concentrations. However, E^c 
does not change. This is in agreement with the published results [26], in which E„c has 
been measured using time of flight measurements. 

4. For a-Si:H(Li), we And that Eo„ increases with the concentration of Li. However, 
Eoc does not change, being 30±5meV for undoped as well as doped samples. This is in 
qualitative agreement with our P doped samples and also with the literature [1,3]. The 
density of mid-gap states is found to increase with Li concentration, again in qualitative 
agreement with our P doped samples (see Table 3.1). This leads us to conclude that 
Li and P are qualitatively similar dopants as far as the bulk states are concerned. 

5. We have obtained the position of D“ in a-Si:H(Li) by combining sub-gap absorption 
and activation energy data. We And that our data favour Street's model of DOS ’9] 
for a-Si:H(Li), in which the position of D" in n type a-Si:H is same as that of D‘^ in 
undoped a-Si:H. This is explained as below. In the other model (by Kocka et al), D" 
moves deeper m the band gap because of intimate pair formation between dopants and 
dangling bond states. However, Li is interstitial dopant and may not form the intimate 
pair. Hence, position of D" in a-Si:H(Li) is the same as that in tmdoped a-Si:H. 
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Chapter 4 


Surface States 


In Chapter 3, we described the determination of the bulk states in undoped and doped 
a-Si:H from sub-gap absorption measurements. This Chapter presents the determin- 
ation of the surface states using these meaisurements. In the literature [1-3], a{hv) 
measurements have been used to determine the density of surface states by measuring 
for a number of samples with identical bulk states but varying thicknesses. 
This is difficult to achieve and gives large uncertainty in the measured value of density 
of surface states. This difficult^' can be overcome by using both CPM and PDS tech- 
niques on the same sample, since the two techniques have complementary informations 
about the bulk and surface states. Smith et al [4] have calculated surface states using 
the relation; Integrated a attributed to the surface states = (Integrated Qr)pi 5 - 2 x 
(Integrated a)c.pm- The factor 2 comes from the assumption that CPM is sensitive to 
about the half the bulk state transitions compared to those seen by PDS. However, this 
assumption is found to hold good only for good quality a-Si:H [5]. Also, this method 
can not be applied to the doped samples. Amato et al [6] use the difference in the 
height of the interference fringes in a(hi/) measured by CPM and PDS to determine 
the density of surface states. However, this method needs perfect uniformity of the film 
and requires identical areas of the film to be probed by PDS and CPM. This is difficult 
to achieve. 

Equivalent occupied surface state density for undoped and doped a-Si:H have been 
measured by Winer et al [7], using total-yield photoelectron spectroscopy. They find 
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that the density of surface states depends on the doping level as well as dopant. For 
B doping, the density of surface states decreases with the increasing doping concen- 
tration. On the other hand, on P doping the density of surface states increases as the 
1/4 power of the doping concentration. Recently Fejfar et al [8] have su.ggested abso 
lute constant photocurrent measurement for estimating the difference in capture cross 
sections between the bulk and the surface states, but do not measure them densities. 
In this Chapter, we suggest an alternative method of measuring surface states, winch 

requires only one sample. 

We have used sub-gap a measured by CPM and PDS to determine the density of 
surface states, a measured by PDS is larger than that measured by CPM in the sub-gap 
region. This difference arises because CPM is not sensitive to the transitions involving 
some of the bulk states and all of the surface and the interface states [4], whereas, PDS 
is sensitive to all the transitions. By measuring PDS for front and rear illumination, at 
high modulation frequency; so that the thermal diffusion length is ot the ordei oi the 
film thickness, Chahed et al [10] concluded that the additional sub gap absorption in 
PDS below leV comes essentially from the surface region. 

In this Chapter, results on the density of surface states for undoped a-Si:H and a- 
Si:H doped with various concentrations of phosphorous and lithium are presented. The 
motivation for studying Li doped a-Si:H is two fold. First Li is interesting because it is 
an interstitial dopant in a-Si:H [11], unlike P which is substitutional. .^Iso. Li has been 
doped by in diffusion of Li in nndoped a-Si:H, at low temperatures (100-I50K). Thus 
we expect that the hydrogen concentration is unchanged upon Li doping. This allows us 
to study the effect of dopant alone. Further, although the bulk stares in a-Si:H(Li) have 
been studied [12], there is no report on the measurement of surface states in a- Si:H(Li), 
so far. 

The present method of calculating the density of surface stares hais the following 
advantages over the existing methods [4,5]. This method requires only one sample and 
there is no specific requirement of the uniformity of the films. Since the technique is 
based on detailed calculations, the assumption that CPM sees about half the bulk states 
compared to PDS is not required. Finally, the technique has been applied to undoped 
as well as n type doped a-Si:H. For CPM, the only requirement is that the sample 
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should be photoconducting. 


4.1 Methodology for Obtaining Surface States 

4.1.1 Undoped a-Si:H 

Chapter 3 gives the methodology to obtain distribution of bulk DOS in the entire 
mobility gap of a-Si:H. Now, in the mid gap absorption region (region III, Fig. 3.3). 
transitions T 4 (Fig. 3.1) are seen by PDS but ignored by CPM. Since, we have already 
obtained the bulk DOS parameters, cx[T 4 ) can be calculated using equation similar to 
Eq. (3.15), from Gaussian distribution of mid gap states. From Eqs. (3.6) and (3.7). 
we get: 

a{Ts) = Opds - ttcpm - Oi{T 4 ). (4.1 ) 

As the distribution of surface states is not known, we have chosen several trial distribu- 
tions of surface states. These include a delta function, a rectangular distribution and a 
Gaussian distribution. All these are assumed to be located near the center of the band 
gap. This assumption is justified because a{Ts) is non zero only in the region III (Fig. 
3.3). Further, we have assumed that the surface states are located below the Fermi 
level. This is done for the simplicity in calculation. Our calculations show that even if 
the surface states are situated differently with respect to the Fermi level, only a small 
difference (<10%) in the value of Ns is observed, provided they are located around the 
middle of the gap. 

From a{Ts) the density of surface states [As(c^~^)] is calculated as follows: 

a{Ts, hu) = {Cjhv) J gi{hu - t)gj{e)de, (4.2) 

where g^ and gj are the densities of initial and final states respectively. The constani 
C includes the transition matrix elements, these are assumed to be the same for all the 
transitions considered here. 

We take 
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gf{i) = 6.7 X 10^' 


(4.3) 

where, e is the energy of the final states measured from the i>an(i rvlges 

The following distributions of surface states, (ji, are used; 

i) 5 function: gi{hi/ - t) = S{e - {hu - Ej)){N,)l{d) 

ii) Gaussian distribution: gi{hiy - e) = {NGld)exp{-ic + E' - hur!'lw- ) 

iii) rectangular distribution: gi{hu - e) = (Nr/d) 

where, Noicm is the height, w the half width at half maximum and E' the 

energy position of the peak of the Gaussian distribution. A;f cm--rr-» ) i.s the height 
of the rectangular distribution and d is the film thickness in cm. 

For the delta function shape of the density of surface states, using Ecps. (4.2) and |4.3) 
we get 


7V,(cm-') = 0(7),) hu d/{K (/m - 
and for the rectangular shape, we get 


1.4) 


Nr(cm-hV-^) = 2a(7;) /m d/[ 3 /v , 4 . 5 ) 

The constant K depends upon the tnatrix element of transition and is assumed to 
Be J. to that obtained for VB extended states to CB extended states transitions 

eooirrr^ir ^ * ^ 

ad St between the calculated o(r,) and experimentally obtained o, r. | i ,| is 

obtatned by varying No. E' and w. This gives .V.. ■ ' . I ■ ■ I 

4.1.2 n-type doped a-Si:H 

Section 3.2.1 describes the determination of bulk DOq ’ fK 

a-Si:H. In the region II (Fig 3 ) „e d fi ' 

(Rg. 3.), we define o,, and by the relations: 


(extrapolated value of a j frnm • tn 

01 apd, trom region I) and 

we write in an obvious notation: 


from region I). Hence ' of 




“pds — <^(^ 4 ) + Oc{Ts) 


aiid, = a(T4). 

Therefore, a{Ts) = is attributed to the absorption involving the surface 

states. Density of surface states can be obtained from a{Ts) using an analysis similar 
to that for imdoped a-Si:H described in Section 4.1.1. 

4.2 Results and Discussion 

4.2.1 Undoped a-Si:H 

Fig. 4.1 shows Oi(Ts) plotted as a function of hv obtained experimentally along with 
those obtained theoretically using a delta function, rectangular and Gaussian distribu- 
tions of the surface states. The analysis is given in section 4.1.1. We find; 

i) For delta function distribution of the surface states: N* = (4.0 ± 1.0) x 10^^cm~^ 

ii) For rectangular distribution of the surface states: Nr is found to be « 6.0 x 10^^ 
{cm~‘^eV~^) and b « 0.6eV; this gives N, = (3.6 ± 1.0) x lO^^cm"^. 

iii) For Gaussian distribution of the surface states: Ng is found to be ft; 1.0 x 10^^ 
(cm~^eV~^)^ w = 0.16eV and E' = l.OeV. This gives = (4.0 ± 1.0) x lO^^cml^. 

These Nj values, for the three distributions of the surface states, are quite close and 
are also in agreement with the values reported in the literature (1 — 5 x 10^^ cm"^) [2,31. 

4.2.2 Li doped a-Si:H 

Using a{hi/) data shown in Figs. 3.11 and 3.12, the density of surface states {Ns) is 
calculated following the analysis of undoped a-Si:H (Section 4.1.2). Calculated Q:(rs ; 
using the three shapes of the surface states along with experimental oi{Ts) are shown in 
Fig. 4.2 for the two Li doped samples studied. We find that the fit is reasonable in all the 
cases. The fit parameters are given in Table 4. 1 . iNL is found to be (S.Oil.O) x lO^^cm 
for sample #1 and (5.5 ± 1.0) x lO^^cm"^ for sample #2 for all the shapes considered. 
Therefore, the density of surface states obtained does not appear to be very sensitive 
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Figure 4.1: a{T.) plotted ae a functiou of h. for undoped a-Si-H(-l SI , , 

"•ith the theoretical lit for Gauasian ffull V ^ ^ ^ 

and r! If f • / (dot dashed line, - . -) 

and delta function (dashed lin^ ^ j- i. n • 

3 ^ . surface states. The error bars 

shown m this figure are the same for Figs. 4.2 and 4.3. 



to the shape of their distribution. These numbers also compare favourably with the 
reported Ns values for a-Si:H(P) [2]. 

So far, we have considered the ’’Standard Defect Model” [15] for bulk states, which 
assumes a narrow Gaussian peak near the mid gap. These are charged defects which 
are all below the Fermi level, in n type material. We now consider the ’’Defect Pool 
Model” [16], in which there may be unoccupied defect states above the Fermi level. The 
optical transitions from the valance band extended states to these states are detected bv 
PDS but ignored by CPM. However, the number of these states are about two orders of 
magnitude smaller than the charged D~ defect states [16]. An analysis, similar to that 
done by us in the case of undoped a-Si:H [17] (Section 4.1.1), shows that the densitt' 
of surface states calculated using the defect pool model is lower by about 2%, which is 
well within the error bars. 

Although we have tried to make the Li distribution uniform by extended annealing, 
a non uniform distribution of Li can not be ruled out, because SIMS profiling was not 
done. Street and Winer [18], based on their SIMS mecisurements, find that surface 
concentration is more than bulk concentration. Since, our analysis assu me s a uniform 
distribution of Li it is worthwhile discussing the effect that the non uniformity of Li may 
have on our results. A non uniform Li distribution may change the density of surface 
states. Also, the DOS in the bulk may vary as a function of depth. Since, our analysis 
assumes the interior states to be uniform, the interior DOS obtained from analysis 
represents an average. However, the calculated surface states density will still be, the 
actual value for that sample, because of the way we have treated the two absorptions. 

4.2.3 P doped a-Si:H(P) 

Fig. 4.3 shows the calculated a(Td for the three distributions (delta function, Gaussian 
and rectangular distributions) of the surface states along with the experimental cx(Ts )- 
The fits are reasonable in all the cases. The fit parameters and the calculated values 
of Ns are given in Table 4.1. We find that the density of surface states under the 
assumption of delta function distribution of the surface states has a value (4.8± 1.0) x 
lO^^cm-^ for sample #5 (E^ = 0.30eV) and (4.2± 1.0)x lO'^'cm-^ for the sample #15 
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Table 4.1: Surface state parameters and density of surface states calculated assum- 
ing delta function, Gaussian and rectangular distributions of the surface states for 
a-Si.H(Li) with two Li concentrations and a-Si:H(P) with two P concentrations. 


Sample 

Details 

Delta 

function 

Ns 

xl0~^^ 

(cm~^) 

Gaussian 

Rectangular 

w 

(eV) 

Ng 

X 10-13 

E' 

eV 

N. 

X 10-13 

b 

(eV) 

N. 

XlO-13 

(cm-3el/-i) 

N, 

xlO-13 

(cm-3) 

a-Si:H(Li)(#l) 

5±1 

0.14 

4.4 

1.0 

5±1 

0.6 

8.0 

4.8±1 

a-Si:H(Li)(#2) 

5.5±1 

0.14 

4.9 

1.0 

5.5±1 

0.6 

9.0 

5.4±1 

a-Si:H(P)(#5) 

4.8±1 

0.14 

4.4 

1.0 

5.0±1 

0.6 

8.0 

4.8±1 

a-Si:H(P)(#15) 

4.2±1 

0.14 

4.0 

1.0 

4.4±1 

0.6 

7.0 

4.2±1 


I 

I 

(Eg. = 0.28eV), independent of the shape of the surface states. I 

I 

As found for the undoped and Li doped a-Si:H, we see that for a-Si:H(P) also. | 


the technique is not sensitive enough to distinguish between the various shapes of the 
surface states. The values of Nj for a-Si:H(P) is found to compare favourably with 
that given in literature [2] (1-5 xlO^'cm"^). These values are also similar to those 
obtained for our a-Si:H(Li). This leads us to believe that a-Si:H(Li) and a-Si:H(P) 
behave similarly as far as the bulk and the surface density of states are concerned. 
We also conclude that the density of surfaces states for doped samples may be slightly 
(within error bars) larger than that in undoped a-Si:H and it does not depend much on 
the doping concentrations. 

4.3 Summary and Conclusions 

Using (Xcpm and ap^s measurements, we have obtained the density of surface states for 
a-Si:H, a-Si:H(P) and a-Si:H(Li) with \-arious concentrations of P and Li. This method 
of calculating Ns has several advantages over the method proposed in the literature, 
such as, this metliod needs only one sample, can be applied to both undoped and doped 
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Photon energy (eV) 


^ "-SPH samples, #5 ( 

nfsl e '■“'angular (dashed line. ---) distril 

oft.sesnrfaces.ates. ^>0- note the different scales for the two samples 




samples and has no specific requirement like perfect uniformity of the film . It, however, 
requires that the sample should be photoconducting, so that CPM can be done. 

N, is found to be about (4.0 ± 1.0) x for undoped a-Si:H. This value of N, 

is in agreement with the literature [2,3]. Further, we have studied the effect of shape of 
the distribution of the surface states on N,. We have taken three distributions of surface 
states namely, a delta function, a Gaussian and a rectangular distributions. We find 
that values of are within error bars of one another, for all the three distributions. 
This shows that the density of surface states is not very sensitive to their distribution. 

The proposed technique for determining density of surface states has been applied 
to Li doped a-Si:H. We find that in a-Si:H(Li) is about (5.0il.0)10^^cm''^ and does 
not seem to vary much for the two Li concentrations studied. Also, is found not to 
depend very much on the distribution of the surface states. However, N* is found to 
increase slightly from the value for the undoped a-Si:H. N* values for a-Si:H(Li) are in 
agreement with the value reported in hterature for a-Si:H(P) [2,3]. 

The density of surface states for a-Si:H(P) is found to be (4.5±1) lO^^cm"^ inde- 
pendent of the doping concentrations for the two concentrations studied. The values 
are in agreement with our values for Li doped samples and also in agreement with 
the literature for P doped samples. Therefore, we conclude that Li and P give the same 
surface states in spite of Li being an interstitial dopant unlike P which is a substitutional 
dopant. 

Finally, we conclude that for all the samples, Nj is not very sensitive to the shape 
of the distribution of the surface states. increases shghtly on Li doping but does 
not depend on the Li concentration. Also Li and P dopings seem to affect the densit\^ 
of surface states in a similar manner. 
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Chapter 5 


Effect of Surface Treatment 


The characterisation of the surface in a-Si:H is important for both understanding of the 
properties of this material and for device fabrication. In particular, it is desirable to 
find a surface treatment steps which reduce the density of surface states. In crystalline 
silicon, this is achieved by cleaving the Si surface in ultra high vacuum, or by heating or 
sputtering followed by annealing [1]. None of these methods can be used for a-Si:H. This 
is because heating at high temperatures (more than 300°C) leads to hydrogen evolution 
from the material and increases the bulk density of states in a-Si:H, quite significantly 
[2]. This renders the material useless, for device purposes. .A.lso, sputtering increases 
the electronic states significantly [3]. Therefore, it is desirable to devise an optimised 
surface treatment step to produce a clean surface having lower density of surface states. 
Attempts have been made to characterise the surfaces for contamination in order to 
produce atomically clean surface in crystalline silicon [1] as well as amorphous sihcon 
[4]. However, no systematic studies of estimation of density of surface states have been 
reported after such treatments, in the case of a-Si:H. 

In this Chapter, we use the technique, described in Chapter 4, for the measurement 
of density of surface states in undoped a-Si:H to see how the density of surface states 
changes when it is subjected to various surface treatments. This also checks the sens- 
itivity of the technique proposed in Chapter 4. We have also calculated the density of 
bulk states (Nj), after various surface treatments using the following relation [5]: 
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iVfe = 1.6 X 10^® / ociifjdE, (5.1) 

J regionll 

wliGrc, ccdt f f tiiG diffcrGiiCG bctwcGii ni63/Sur6(i Of using CPM, in the region II 3 ,nd 

extrapolated value of a from region I (Fig. 5.1). 

5.1 Surface Treatments 

We have used various surface treatment steps, which are defined below. .State .A. is the 
rested state of the as deposited sample. States B and C are the light soaked ( 20 m W/ cm^ 
for 15 hours) and the a nn ealed {15Q°C for 2 hours) states respectively. 

State D is obtained by etching (HF etched. 15 % in DI water) the native oxide from 
a-Si:H surface. The surface is again oxidised by treating with concentrated H'^SO^ 
and subsequent annealing (State E). States F and H are the room temperature plasma 
treated states. State F is obtained by treating the sample in state E with room temper- 
ature nitrogen plasma (300K, AQmSN jew? for 2 hours). State G is obtained by annealing 
the sample in state F (Annealed at 150°C for 2 hours). State H is hydrogen plasma 
(300K, 40mW/cm- for 2 hours) treated state. More details of the various treatment 
steps are given in Section 2 . 1 . 4 . 

5.2 Results and Discussion 

In Figs. 5.1 - 5.4, we present our results on sub gap absorption measured as a function 
of hi/ by CPM and PDS for an undoped a-Si:H (Sample No. #18) in state A (rested 
state) and state B (light soaked state) [Fig. 5.1]; States D (oxide etched state) and E 
(oxidised state) [Fig. 5.2]; States F (nitrogen plasma treated state) and G (subsequently 
annealed state) [Fig. 5.3]; and state H (hydrogen plasma treated state) [Fig. 5 . 4 ]. Sub 
gap absorption data is analysed and the density of surface states (N,) are obtained 
using the analysis described in Chapter 4 [6,7j. Also, the density of the bulk states (Nj) 
are calculated using Eq. (5.1). The values of and N 5 are plotted in Fig. 5.5 after 
various surface treatments; from states A to H. 
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hi;’ ( e V ) 

Figure 5.1: log a as a function of hz/ as measured by CPM, PDS and Transmission 
for a-Si:H sample. Fringes in the curves arising from the thin film interference have 
been averaged out. States A and B represent the rested and the light soaked states 
respectively. Note different a axes for the two states in Figs. 5.1 to 5.3. The error bars 
in this figure is the same for the Figs. 5.2 and 5.3. 
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Fig. 5.1 shows a(hz/) measured by CPM and PDS in the rested state (state A) and 
after light soaking (state B) in undoped a-Si:H. and Nj calculated using CPM and 
PDS data in Fig. 5.1, are shown in Fig. 5.5. From Fig. 5.5, we find that after light 
soaking (State B) the density of dangling bond states (Nj) increases from 5 x lO'^-m-^ 
to ~ 1 X 10^^cm“^ without much change in the surface states. This is in agreement 
with the literature [8,9]. Annealing the sample (150“C for 2 hours; state C) brings 
back to its rested value (« 5 x lO^^cm"^). 

Fig. 5.2 shows a{hi/) measured using CPM and PDS on etching a-Si:H in state C 
(annealed state) with 15% HF in deionised water (state D). Fig. 5.2 also gives a(hu) 
measured using CPM and PDS after re growth of oxide layer (state E) by treating the 
sample with concentrated 112504 followed by annealing. We find (Fig. 5.5) that the 
density of surface states increases from ~ 4xl0^^cm~^ in the annealed state (state C) 
to « 6.5xl0^^cm"^ after etching with HF (state D). The bulk states, however, remain 
unchanged within the errors. This shows that freshly etched surface, which is de\'oid 
of the native oxide, has a higher density of surface states. This may be explained as 
follows. One type of structural disorder that might lead to the surface states in the gap 
is the presence of near surface dangling bonds [10]. Removal of oxide from the surface 
might increase the near surface dangling bond and hence an increase in the density of 
surface states on oxide etching. We confirm this by growing oxide again and measuring 
the surface states. The oxide passivates the surface and brings back the density' of 
surface states to its starting (lower) value in the annealed state (state E). 

Finally, we have treated the a-Si:H sample first with nitrogen plasma (state F) and 
then with hydrogen plasma. Fig. 5.3 shows a{hu) measured using CPM and PDS in 
state F along with that in state G. State G is obtained by annealing the state F. Nitrogen 
plasma (300K, AOrnW / cm^ ,dc plasma for 2 hours) treatment is found to increase the 
density of surface states from ~ 4 x lO^^cm to 5 x 10^^cm~^. This small change 
is within the errors of measurement. Also there is a small increa.se in the density of 
bulk states, which is within the errors of measurement. This shows that treatment 
with nitrogen plasma does not have much effect on either the bulk states or the surface 
states. No change in bulk states as well as the surface states was observed on annealing 
the sample after nitrogen plasma treatment. 
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hV (eV) 


Figure 5.2: log o; as a function of hi/ as measured by CPM, PDS for a-Si:H sample. 
Fringes in the curves arising from the thin film interference have been averaged out. 
States D and E represent the HF etched and the oxidised states respectively. 


on 
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Figure 5.3: log a as a function of hi/ as measured by CPM, PDS for a-Si:H sample. 

Fringes in the curves arising from the thin Hm interference have been averaged cut. 

States F and G represent the nitrogen plasma treated and the subsequent annealed 
states respectively. 




Fig. 5.4 shows a{hv) measured using CPM and PDS after hydrogen plasma treat- 
ment (300K, 40mW/cm^dc plasma for 2 hours). Hydrogen plasma is found to do the 
surface passivation by decreasing the surface states to « 2.5 x 10‘^cm"^. This may 
be because of surface dangling bond termination by H atoms in the plasma. We thus 
conclude that Hydrogen plasma gives a surface with a low density of surface states 

(» 2.5 X 10‘^cm-^). Etching of native oxide, however, increases the density of surface 
states by about 50%. 

5-3 Summary and Conclusions 

We have measured sub gap absorption coefficient as a function of photon energy us- 
ing CPM and PDS for a-Si;H after various surface treatments. These include light 
soaking (B), annealing (C), etching of the native oxide using HP (D) and subsequent 
oxide growth by concentrated sulfuric acid (E) and treatment with nitrogen (F) and 
hydrogen plasma (H). Density of surface states and the bulk states are calculated after 
each surface treatment. We find that light soaking increases the density of bulk states 
by a factor of 2, without changing the surface states. This is in agreement with the 
literature that light soaking is mainly a bulk effect. Another notable results regarding 
the surface states are that the surface states increase from w 4xl0^^cm~^ in the an- 
nealed state (state C) to ~ 6.5xl0^^cm“^ after etching with HF (state D). Hydrogen 
plasma treatment, however, is found to decrease the density of surface states from ~ 
5x 10^'^ cm“^ in state G) to ~ 2.5xl0^^cm~^ after etching (state H). This explains why 
exposing with hydrogen plasma produce better results whenever a-Si;H based devices 
involving several interfaces (e.g. solar cells) are made. 
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Figure 0.4: log a as a fimction of in as measured by CPM, PDS for a-Si:H sample. 
Fnoges in the curves arising from the thin film interference have been averaged out. 
States H represents the hydrogen plasma treated state 
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Figure 5.5: Bulk and surface density of states after various surface treatments. State A 
(Rested), State B (Light Soaked, 2 {imWlcw? for 15 hours). State C (Annealed at 150°C 
for 2 hours). State D (HF etched, 15 % in DI water). State E (concentrated. H 2 SO 4 : 
treated and annealed). State F (exposed to Nitrogen plasma, 300K, 40mW/cm^ for 2 
hours). State G (Annealed at 150“C' for 2 hours) State H (exposed to Hydrogen, 300K. 
40mW/cm2 for 2 hours). The vertical lines in the state A gives the error bars of the 
measurements (after reference 6). 
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Chapter 6 

Potential Fluctuations 

6.1 Potential fluctuations in a-Si:H 

So far, we have explained sub-gap absorption results in terms of single particle densit\' 
of states, assuming a-Si:H to be a homogeneous material. In reality, however, a-Si:H 
is heterogeneous. The heterogeneities give rise to potential fluctuations in a-Si:H. Sev- 
eral experiments like NMR [1] show that even the best quality (device grade) a-Si:H 
contains voids and clusters of hydrogen, giving rise to potential fluctuations. Addition- 
ally, potential fluctuations may also arise from the variations in bond length and bond 
angles. Further, neutral and charged dcingling bonds as well as the dopants (if present ) 
may also contribute to the potential fluctuations in a-Si:H [2,3]. It is not easy to deal 
with the potential fluctuations in a fundamental manner. Nevertheless, it is important 
to study them because they are expected to influence the electronic properties in a 
profound manner [4]. 

The range of these potential fluctuations has been estimated to be of the order of 
inter atomic distances (~ 3A) [5]. However, in the transport measurements (conduct- 
ivity and thermopower) the carriers move through these potential fluctuations. So the 
carriers are likely to overlook very sharp changes in the potential because of quantum 
mechanical tunneling. Therefore, these measurements are sensitive to the range of fluc- 
tuations, which is likely to be longer (^ 80A) [6,7]. Overhof and Beyer [6,7], have ex- 
plained the difference between the acti^-ation energies of conductivity and thermopower 


95 



by coasidering the long range potential fluctuations [6,7|. The details are described 
below; 

The general expression for conductivity, assuming extended state conduction by 
electrons, is given as 

cr = j a{E)dE = aoexp{-{Ec - Ef)/kT) = a„exp{-EJkT). ( 6 . 1 ) 

where, a, contains the mobility factor. E, = E, - E; is the conductivity activation 
energy. Similarly, thermopower S can be written as 

S = -^j{E-EMB)dE, ( 6 . 2 ) 

For the extended band conduction, Eq. (6.2) gives 

S = --[(£. - E,)lkT + .4] = --{BJkT + .4). (6.:!) 

e € 

where, A 1, is a constant. Again, E, = Ec - E/ is the thermopower activation tmergy. 
In a-Si;H, in general we expect = Ej, contrary to the observations. The following 
reasons have been proposed for this dilference [8]: 

1. If both electrons and holes are conducting. S will have opposite signs for the two 
carriers, making Ej different from E^y. But this effect can be neglected for doped 'n- 
type) a-Si:H, where, only one type of the carriers (electrons) dominate. 

2. Mobility edge may not be very sharp. 

•3. Conductivity may not be spatially homogeneous, because of the heterogeneities in 
the material. This appears most likely in a-Si:H. 

In order to eliminate the effect of Fermi level position because of either the statistical 
shift or any other agency like doping, Overhof and Beyer [4,5], have defined a function 
Q given by 

Q=lna + {elk)S = Qo + EQ/kT, (6.4) 

where, =ln tJo + A and Eg = E^r - E^. In order to explain non zero value of Eg, 
Overhof and Beyer [6,7] have done model calculations for c and 5, in which a-Si;H 
sample is divided into small cells. Each cell has a disorder potential attributed to it by 
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random number generation, which decides the local mobility edge. Their model shows 

that Eg is a measure of the width of the long range potential Unctnations (A), and the 
two is related by [6,7] 

Eq = 1.25A (6.5) 

As already explained, these measurements can give only the long range potential 
fluctuations, which may be some kind of average of the short range potential fluctuations. 
Optical measurements, on the other hand are the vertical transitions. Therefore, optical 
measurements are likely to sense the short range potential fluctuations, inherent in the 
material. Attempts [5,9] have been made to calculate o.ijiu'j in the sub gap region 
using potential fluctuations, theoretically. These attempts are limited to obtaining the 
exponential absorption in the Urbach edge region, starting from the Gaussian potential 
fluctuations. It is found [10] that for the range of potential fluctuations, L » A (where 
A is the de broglie wave length of the free electrons); a Gaussian Urbach edge is obtained 
from a Gaussian potential fluctuations. On the other hand, Halprin and Lax [11] find 
that for L << A, the density of localised states in the Urbach edge, g(E) oc exp{E^^'^). 
Quantum mechanical calculations of John et al [5], successfully predict the exponential 
shape of the Urbach edge, starting from the Gaussian potential fluctuations, provided 
the correlation length of the potential fluctuations is taken to be ss 3 .4. Thus, we expect 
that the range of the potential fluctuations seen by optical absorption measurements are 
short, of the order of the inter atomic distances (~ 3 A). The calculations by John et 
al [5] show that the width of the short range potential fluctuations (Vj) is related to the 
slope of the Urbach edge £„ by the relation [5]: 

V; = 20EoeL, ( 6 - 6 ) 

where, ci = j L'^) ~ 0.5eV. 

6.2 Determination of potential fluctuations 

In this section, we first determine the width of long range potential fluctuations using 
transport measurements in undoped a-Si:H and a-Si:H doped with various concentra- 
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tions of phosphorous (P). We then estimate the short range potential fluctuations using 
sub gap absorption data, measured by PDS on the same samples for which width of 
long range potential fluctuations were estimated. An empirical relation between the 
two widths is discussed. Also, we explain the difference in o(Jn>) mea.sured by CPM 
and PDS qualitatively. 

6.2.1 Determination of potential fluctuations using transport 
measurements 

We estimate the long range potential fluctuations from the transport data following 
Overhof and Beyer [6,7]. We take to be constant in a cube of linear dimensions L (% 

80 A). It IS, therefore, possible to define a local conductance inside a cube of volume 
centered around a site X,- as: 


(cTal L)exp[ {Ec(Xi) — Ep')fkT] 

e, Ec(W,) - < Ec > +V(X,-) IS the position of the local conduction band. The 
probability that V(Xi) has a value e is given by [6,7] 






Therefore, the sample is replaced by a resistance network, and gives 0 = i„teroo„- 
oectmg resistances, , ,,, ^ 

mo el calculations. We use Kirchoffs relations to solve the problem, 
e, nations are solved by theinversion of „3 , „3 

the hL^rT"; T*’" 

of 1/T pl T H '= S) as a function 

ot 1/ 1 plot IS related to the width X as 


Using Eq. (6.9) 
Q function plotted 


we calculate the width of potential flnctnations (A). Fig. 6.1 shows 

“ ^ of >000/T for three P doped samples. E, is estimated 
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Table 6.1; Conductivity activation energy (Eo-), slope of Q function as a function of 1/T 
(Eg) and width of the long range potential fluctuations (A). Samples are listed in the 
order of increasing P concentration. *Eq for undoped sample is taken from Ref 12 and 
is not for the sample #18 


Sample 

details 

E. 

(eV) 

Eq 

(eV) 

A 

(eV) 

a-Si:H(#18) 

0.70±0.01 

(.05 ±.01)* 

.04 ± .005 

a-Si:H(P)(#5) 

0.30±0.01 

.08 ±.01 

.06 ± .005 

a-Si:H(P)(#12) 

0.28A0.01 

.08 ±.01 

.06 ± .005 

a-Si:H(P)(#ll) 

0.24A0.01 

.10 ±.01 

.08 ± .005 


from the slope of these curves. A for all the samples is calculated. Values of Eg along 
with A are given in Table 6.1. It may be noted that Q for undoped sample could not be 
measured because the sample is highly resistive. Eg in Table 6.1 has been taken from 
the literature [12] for a similar sample. We find that A increases from 0.040± O.OOSeV 
for the undoped sample (sample No. 18, E, « 0.7eV) to 0.08± O.OleV for the sample 
with highest doping (sample No. 11, « 0.24eV). 

6.2.2 Determination of potential fluctuations from optical ab- 
sorption data 

Fig. 6.2 shows sub gap a{hu) measured between 0.9eV and 1.8eV using CPM [13,14] 
and PDS [15,16]. It is found that in the sub gap region, Op* > a^rn- This difference is 
qualitatively explained using potential fluctuations as follows. 

Fig. 6.3 shows schematic of the potential fluctuations as might be present m a-Si:H 
[17], E, and E. represent the conduction band and the valance band percolation levels 
respectively. Representative optical transitions are also shown m the figure. Transitions 
of type 1 takes an electron to the final state above the conduction band percolation level: 
whereas, transitions of type 2 and 3 take the electrons to the levels below the critical 
percolation level. In CPM. only those transitions are seen, which are able to conduct 
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a-Si: H(P) 

Sample 

#5 


EQ(eV) 

^ 0.08 


d 6 


#12 


0.08 



Figure 6.1: Function Q = Incr +(e/k)S plotted as a function of 1000/T for three P 
doped samples; sample #5, sample #12(ref. 8) and sample #ll{ref. 8). Note the 
different scales for different samples. 







Figure 6.2: Sub gap absorption measured as a function of photon energy using 
and PDS(regions II and III) and using transmission(region I) for undoped a-Si:H. The 
interference fringes in the data have been averaged out. 
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i.e. which are above level. So, CPM is not sensitive fo traissPint,- r , 

^ Ml lypp 2 and 

3. On the other hand, PDS is a thermal effect and is likeiv hv 

tr) all the 

transitions from 1 to 3. This qualitatively explains iln* dilffre- rv 

I and 

acpra m a-Si:H. As estimation of in the potential 

' rMf|i]]pes ^ 

knowledge of both the optical absorption as well as frans,,.,,. i„ 

and is not discussed here. We used PDS data for „.a„„.a„v,. 

fluctuations. ' 

In the literature, considerable amount of work [.u.dl I, as done lo explain th 

exponentialnatureoftheoptical absorption edge in .he I-rl, a, -I, reuu,„. s„ far ,l„.o„r 

successful attempt in explaining the exponential edge starting fro,,, , 

iuctuatrons, seems to be by John et al (5]. Th. autbor.s show ti.a, ,f [ ~ A i, ' 

poss.ble to get an exponential absorption edge althongh p„,e.„,u,| 

a G.ss,an dmtrrbnt.on. Indeed they argue that L . A is the exper.n.entaliv observed 
™ e m tbe Urbach edge regiou. Here, J h,Jx, , 

H:ba°chtr “ in 

uroach edge region is given by 


'’”<'^=A“[-7~eap(-r/vS5 


The slope of Urbach edge E, is related to V 


i as [o , 


{ 6 . 10 ) 


, ^ ( 6 . 11 ) 
where, ei = % O.oeV fSl Pm fi 

using Eq. (6.r ) for all th "^«a«ured Urbach .slop,. y,, i.s ralrulated 

L ^■^ *1“"' samples, and are given i„ Table tU 

i-et US now consider region II rpio- fi 4 W i 
been made in the Uteratnre to calcuLe poltl, 

“ '^Sion- The potential Huctnations in thrs re ' ^ 

dangling bonds. We nronn. fK u 

distribution of potential iuctuations. Theltltl! IT'*' ^ 

*e potential fluctuations given by Eq (6,im w superimposed on 

of the potential fluctuation <; h ' ^ this second component 

nations has a smaller heiah^ . j , 

g and a larger width than the first 
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Figure 6.3: Long range potential fluctuations as might be present m a-Si:H(Ev and Ec 
represent the \^lance and conduction band mobility edges; Ep is the Fermi level). CPM 
sees the optical transitions marked 1, ivhereas, PDS sees all the transitions marked 1 

and 2 (after reference 17). 
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component (V^ » 14). We further assume that the width of the sr. fUKi r,,ni[>on»*iit of 
the potential fluctuations is also given by a relation .similar to that hi . fi.ii); 

{V:Y = 2QEm { 6 . 12 ) 

where, is the slope of the region II of the optica! ahsor|>tion ri]r\'<* c Fig. fi..!). 
Therefore, the probability, p(e), that the local conduction band f away from its 
mean value is given by the following sum of two Gaussians. one farii licrause of the 
absorption from the Urbach edge region and from the plateau rfnion: 

^ \ ■; )-) (6.1.3) 

where the zero of the energy is chosen at the mean of the iluct uutiug <-on(iuclion band. 
Pa is the normalisation factor for p(€) and is found to be p„ = y (2}/'i 1 t-B). We deiiiu' 
the total width (Vr) for all the potential fluctuations taken together by: 

Vt = [ (6.14) 

J ~«x> 

Using p(e) from Eq. (9), we get = [{V^f + 4- The scaling 

factor B is calculated by taking the ratio of qi l.SeU) and the extrapolated value of a 
at l.SeV from the plateau region. 

We find that E, increases with the doping level: being 55meV for the undoped 
a-Si:H (E^ Ri 0.7eV) and si llOmeV for highest doping (E, i).2I<AU. The width 

of the potential fluctuations (Vt), estimated from the optical ab.sorption data, is also 
found to increase with the level of doping, being = 5 ; O.SeV for the uinioped sampl(> (No. 
lb) and 1.45eV for the sample (No 11). with the highest P doping. The values of E' , 
d Vy for all the samples are given in Table 6.2. Also, the width of the long 
g potential fluctuations (A) estimated from the transport measurements is found 
to increase with the doping level. A ^ 0.04eV for the undoped sample (No. 18) and « 
O.OSeV for sample with highest P doping (No. 11). 
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Figure 6.4: Sub gap absorption as a function of photon energy measured using PDS for 
nudoped(#18) and P doped a-Si:H(#5 and #11). The interference fringes have been 

averaged out. 
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Table 6.2: SIope(Eo) of Urbach edge. width{V^) of fhf fif.f 
potential fluctuations, SIope(E') of plateau part of 


'<■ 'b* Uiurf 


range 


-.viritht V' . . ,r . , 

• 011,1 term 


of the short range potential fluctuations, B and total widflnX’-' a' •, 


potential fluctuations. The samples are listed in the ortic 
concentration 




range 


oa ' ' •» JJliu'ipinjrous 


Sample 

details 

Eo 

(meV) 

(eV) 

K 

(eV) 

-j " ^ 

I 

{«‘Vi 

a-Si:H(#18) 

55±5 

0.74i.05 

0.56±.02 

’ 2.;]7 r.().5 

a-Si:H(P){# 5 ) 

85±5 

0.92±.05 

0.51 ±.02 I 

! i 

’ 2.26 ± .nr, ' 

a-Si:H(P)(#12) 

100±5 

1.0±.05 

p 

CO 

hS 

2.19 2: .nr 

a-Si:H(P)(#ll) 

110±5 

1.05±.05 

0.45±.02 1 

1 . 12 2 : .nr ! 


B 


\ J. 


:.!)(! i 


r.uni 


0.7.-) r 
1 ■!,).") - 
1.20 r 
1.00 -t 


.0.5 

.05 

.05 

.05 


6.3 Discussion 

We fed that V, > 4 foe all the samples st.,die,l. This is ,„.r,,r,si„p ,ra,.s- 

ort measmemeots ma. a„t see the shasp sh„„ , ,, 

can tunnel through them (,eantam mecha„KaU,„n,eli„ei riuTef , , 

measurements see only the l„ ' ‘"'''"“P”' 

..2 JZZI""' "■ ■■ 

itive to Short range potentialfl t ■ "e seus- 

■~r..tT j' “ “ '■ “ *"■ — 

y. tar. It .s mteresting to see how the width of , I 

scale with the ranee of the t ■, ‘™“ "f '1“‘ Potemial llm-tuatioiis 

fan„e of the potential fluctuations. We olot V r ■ 

6.5, and And that the data es ». P'°‘ Vy as a function of A in Fig. 

'‘•'data can fit a straight line of the form- 


= KA 

our samples. Qualitativelv. 


(6.15) 


Here, K is found to be 20±5 for 

“ potential fluctuations) wr^riforT""' " 

kind of averaging of the fluctuations. 
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The aver^^ging decreases the width of the Buctuations by a factor of ss 20. in our case. 
Further, it is interesting to see that K ss L/L', where, L and L' are the ranges of the 
long and the short range potential fluctuations, which are 80 A and 3 A respectively. 
This shows that the transport measurements are sensitive to a potential fluctuations 
which get averaged by a factor of ss L/L'. This is, however, only an empirical result 
and needs to be further justified theoretically. This has not been done, so far. 


6.4 Summary and Conclusions 

We have proposed explanation of sub gap absorption measurements using potential 
fluctuation model, in undoped as well as P doped samples. This is an alternative 
explanation and is different from the single particle density of states explanation of 
optical absorption, in which a-Si:H is assumed to be homogeneous. We find that apd^ 
is larger than for all the samples. This result is in agreement with the literature. 
In the potential fluctuations model, this difference is explained as follows: 

In CPM, we meaisure a from photoconductivity. In undoped and n type doped a- 
Si:H, at room temperatures and higher only electron conduction is important. This is 
because the hole mobility is about an order of magnitude lower [20]. Also in CPM, the 
photogenerated carriers must be able to conduct in order to contribute to a^. On the 
other hand, PDS measures the absorption coefficient by measuring the heat produced 
when the sample absorbs light. Thus, only the transitions of the type marked 1 (Fig 
6.3) contribute to CPM, but all the transitions marked 1 to 3 contribute to PDS. This 
explains the experimentally observed result {apds > otcpm) qualitatively. However, it is 
difficult to obtain occpm theoretically. Therefore, for calculations of the width of short 
range potential fluctuations, we use only the PDS data, which is easier to handle. 

We have estimated the width of potential fluctuations in undoped a-Si.H and a-Si.H 
doped with vatic us amounts of phosphorous (P ), using optical sub-gap absorption meas- 
urements and transport measurements, viz conductivity and thermopower as functions 
of temperature. Both measurements show that the width of the potential fluctuations 
increases upon increasing the phosphorous doping. We find that the width of the po- 
tential fluctuation measured by optical measurements is larger than that measured by 
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transport measurements (A) by a factor of « 20, for all the samples. This discrepancy 
is explained by noting that m transport measurements, carriers are likely to overlook 
sharp fluctuations in the potential because of quantum mechanical tunneling, as these 
fluctuations are short ranged, of the order of inter atomic distances. Thus, transport 
measurements see long range potential fluctuations whose width may be smaller than 
that of the shot t range fluctuations. On the other hand, optical absorption is a local 
process and is likely to be sensitive to the short range potential fluctuations inherent 
in the material. We find empirically that the averaging decreases the width of the po- 
tential fluctuations as seen by the transport measurement by a factor of ss 20. It is 
also interesting to note that V^/A ss L/L , where L and L’ are the ranges of the long 
range and short range potential fluctuations. Therefore, we conclude that the transport 
measurement see a potential fluctuation which gets averaged by a factor of L/Lh This, 
however, needs to be justified theoretically. 
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Chapter 7 


Summary and Conclusions 


We have setup constant photocurrent measurement (CPM) [1] and photothermal deflec- 
tion spectroscopy (PDS) [2]. Complete automation of the measurement have been done 
to get reliable set of data. Sub-gap absorption (a), for un doped hydrogenated amorph- 
ous silicon (a-Si:H) and a-Si:H doped with various concentrations of phosphorous and 
lithium, has been measured using PDS and CPM. We find [3], in agreement with the 
literature [4], that a measured by PDS is always greater than that measured by CPM 
on the same sample. This is explained by arguing that PDS is a photothermal effect 
and hence is sensitive to all the transitions. On the other hand, CPM sees only those 
transitions which terminate in the conduction band (CB) extended states. Further, 
unlike PDS, CPM is not sensitive to the transitions involving surface states because 
the carriers generated at surface and interface region quickly recombine [5] without 
contributing to the photo current. 

In the present work, a measured by CPM and PDS on a given sample of given 
thickness is analysed to get bulk states in the entire mobility gap as well as surface 
states [3]. This method obviates the problem of making identical samples with different 
thicknesses as required by the other methods of measuring surface states [6-8]. However, 
the shape of distribution of DOS needs to be known for the analysis. Various shapes 
of the mid-gap states like Gaussian [9] and exponential [10] have been suggested. W'e 
find that CPM data can not distinguish between various shapes of the mid-gap states 
[1]. This conclusion was drawn by Payson et al [11] for their PDS data. We show that 


112 



this ambiguity can be resolved using dual beam constant photocurrent measurement 
(DBCPM) [12]. We find that the anomalous hump appearing in the DBCPM spectrum 
at lower chopping frequencies, can be fitted only if we assume the Gaussian distribution 
for the mid gap states and not the exponential. Thus, we conclude that the mid gap 
states are most likely distributed as a Gaussian. Further, we have also calculated the 
correlation energy (U) between electrons occupying the same defect site for undoped 
a-Si:H, using the DBCPM spectrum at a high chopping frequency (333Hz) and DC 
CPM data. We get U = O.lSeV, in agreement with the literature [12]. By plotting 
logicxpds - cccpm) as a function of hu in the Urbach edge region, we get E^c, the slope 
of the conduction band Urbach edge [3], which is not possible using OTp^s or c^cpm data 
alone. 

We have calculated the bulk states in the case of undoped, lithium (Li) doped and 
phosphorous (P) doped a-Si:H, using sub-gap absorption data. Li is an interstitial n 
t\^pe dopant in a-Si;H [13], unlike P, which is substitutional. We find [14] that the slope 
of the valance band Urbach edge (Equ) increases with the increasing P concentration as 
well as Li concentrations. The slope of conduction band Urbach edge (Eoc), however, 
is found to 30±5 meV, independent of P or Li concentrations. Also, the density of 
mid-gap states are found to increase with the increasing P or Li concentrations. Our 
results for P doped samples are in agreement with the literature [15,16]. No results are 
available, in the literature for Li doped samples, at present. Our findings show that P 
and Li dopants give similar results. 

By combining the sub-gap absorption and the activation energy data, we have ob- 
tained the energy position of the D“ in a-Si:H(Li) [17]. Two distinct models for the 
position of dangling bond states in doped a-Si:H have been proposed. According tC' 
Street et al [9], position of D~ in doped a-Si:H is identical with that in imdoped a-Si:H. 
Kocka et al [18], on the other hand, propose that the position of D~ in doped a-Si:H 
moves deeper in band gap, moving away from the CB edge, because of the intimatC' 
bond formation between the dopant and the dangling bond states. We find that even 
after careful analysis of the CPM spectra, the two models are within the experimental 
error. But on using the conductivity data along with the acpm data, we find that our 
data prefer the Street’s model in a-Si:H(Li) [17]. This may be because Li is an intersti- 
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tial dopant and may not form intimate pairs with the dangling bonds for Kocka’s model 
to be valid. 

One of the inputs required for the analysis for determination of surface states is the 
shape of distribution. Since this is not known, we have taken three shapes of the surface 
states namely delta function, Gaussian and rectangular. The density of surface states 
for undoped a-Si:H is found to be (4 ±l)W^cm'-^ for all the three shapes. Therefore, 
we conclude that the density of surface states is not very sensitive to its distribution 

[ 3 ]. 

The surface states in a-Si:H(Li) hcis been measured for the first time and compared 
with a-Si:H(P) for various concentrations of Li and P. We find [14] that the density of 
the surface states in a“Si:H(Li) is ((5 ± 1) lO^^cm""^), independent of the shape of their 
distribution, and is about 20% larger than that in undoped a-Si:H. It is found to be 
independent of the concentration of Li, for the two concentrations used in this study. 
The density of surface states in a-Si:H(P) is found to be about the same as a-Si:H(Li). 
This shows that density of surface states in the two dopings are similar. 

Further, we have measured the density of bulk and surface states, after exposing 
undoped a-Si:H sample to light and various surface treatments. These are etching the 
surface oxide with dilute hydrofluoric acid, growth of oxide layer by treating the sample 
with concentrated sulfuric acid and subsequent annealing and exposure to nitrogen and 
hydrogen plasma. We find that upon light soaking, the density of the bulk states in- 
creases by a factor of 2, without much change in the surface states. This is in agreement 
with the common conviction [19], that light soaking is a bulk effect. Etching of native 
oxide with hydrofluoric acid, increases N5 by a factor of about 1.5. This shows that 
the surface devoid of the native oxide have a higher densit}" of surface states. Treating 
the etched sample with concentrated H 2 SO 4 to grow oxide again and annealing brings 
back the density of surface states to its value in annealed state. 

In addition, we find that nitrogen plasma treatment increases the density of sur- 
face states slightly. This small change is within the errors of measurement. Hydro- 
gen plasma, however, passivates the surface by decreasing the surface states by about 
50%. This may be because surface etching and the Si bond termination by H atoms in 
the plasma. This might explain why the interfacial surfaces are exposed to hydrogen 
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plasma, while making a-Si:H based devices. 

Finally, we have explained sub“gap absorption results using potential fluctuations 
[20]. This is an alternative explanation, different from the single particle density of 
states explanation, in which a-Si:H is assumed to be homogeneous. The total width of 
potential fluctuations (Vt) is estimated by optical absorption measurements, using PDS. 
The two regions (Urbach edge region and mid-gap absorption region) of the optical ab- 
sorption curve have been explained by proposing two Gaussian terms in the potential 
fluctuations. Also, the width of long range potential fluctuations (A) is calculated from 
transport measurements, for the same samples. We find [20] that the width of the poten- 
tial fluctuations obtained from optical as well as from transport measurements increases 
with the increasing phosphorous concentration, as expected. Also, we find that Vy is 
larger than A by a factor of 20, for all the samples. This discrepancy is explained 
by arguing that optical and transport measurements are sensitive to different ranges of 
the potential fluctuations. This is in line with the argument of Fritzsche [21] explaining 
the difference between electrical and optical gaps in amorphous semiconductors. We 
find an empirical relationship (Vt ~ ATA) between the long range and the short range 
potential fluctuations, where, K = 20±o. This shows that the transport measurements 
see potential fluctuations, which are some kind of average of the short range poten- 
tial fluctuations, seen by the optical absorption measurements. This emperical results, 
however, needs to be verified theoretically. 


Scope for Future Work 

This work demonstrates the utility of the sub gap absorption measurements for obtaining 
information about the bulk as well as the surface states. An obvious extension of the 
present work is to use the technique in the case of boron doped a-Si:H. Studies of the 
surface states after various treatments might give the reason for lower density of surface 
states in a-Si:H (B) compared to that in undoped a-SicH. Similar studies on compensated 
a-Si:H should be interesting. Also, effect on density of surface states of alloying a-Si:H 
with other materials like germanium, carbon, tin etc. Some of these materials are useful 
for devices and studies of their electronic properties are desirable. We have explained 
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the difference in and in terms of the optical transitions not seen by CPM 
but seen by PDS. We find that in a chalcogenide glass (SegoTe^o), both and 
are the same [2]. Further experimental work and analysis is necessary to explore and 
explain this observation in chalcogenide glasses. PDS is a unique tool for probing non 
radiative transitions. Luminescence, on the other hand, gives the radiative transitions. 
When these two techniques are applied on the same sample, a more complete picture 
of the branching ratios of radiative and nonradiative transitions in the material imder 
study, may emerge. 

dcpmihu^ have only been explained qualitatively in this work, using potential fluc- 
tuation model. CPM requires knowledge of both optical absorption and transport of 
carriers through the potential fluctuations, for explanation. This requires a better the- 
oretical understanding of these ideas. This also is a possible direction for the future 
research. 

This work is only an initial effort in explaining sub gap absorption using potential' 
fluctuations. Many problems remain to be solved. For example, the relation between the 
long range and short range potential fluctuations, found in this work, is only emperical. 
A theoretical justification of this relation requires the understanding of the averaging 
processes, when we visualise the long range potential fluctuations as an average of the 
short range potential fluctuations. If achieved, this will allow the model to produce 
transport properties, from the optical absorption data. This would be of fundamental 
as well as practical interest. In conclusion we can say that a lot of unsolved problems, 
controversies, issues and ideas remain to be explored in a-Si:H and its alloys. 
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Appendix A 


Transmission Measurement 


Under the assumption of : (a) non absorbing substrate and (b) K=0 which is true for 
most part of the spectrum 


where, 


T = Axl{B — Cx.cos(f) + Dx^), 


(A.l) 


A = 


B = {n + l)^(n + s^) 


(A.2) 

(A.3) 


C = 2(n^ — l)(n^ - s^) 

(A.4) 


{A.5) 

$ = iTrndfX 

(A.6) 

X = exp{—ad) 

(A.7) 


The symbols are defined in the Fig. Al. Fig. A2 Shows the transmission spectrum of 
a-Si;H thin film. Tm and T,„ are envelopes of maxima and minima of the interference 
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fringes. The spectrum is divided into three regions : transparent region (r„ > 0 . 91 ) 
weak and medium absorption region {0.91 > I),, > 0.5) and strong absorption region’ 
{Tm > 0.5). In these regions; n and o are given in the closed form [Ij. 

(a) Transparent region: 

n = + 32)1/211/2 

■' (A. 8) 

where 

where, M= 2s/Tm - + l)/2 

and a = 0 

(b) Weak and Medium Absorption Region: 


n = 


= (A + (/V2 -52)1/2)1/2 
where, N = 2 x + sl±l 


(A.9) 


and X = ^ ~ — 3^)]i/2 

{n - l)3(n - 52 ) 

f 

where P — 

’ nprZ 


(A.IO) 


W ^strong Absorption Region: As seen from fig. A2. there :s no interference fringes 

n(A) -Si™- "( A) is calculated by the extrapolation of 

ri;;,: ^ - ^-nm and transparent) using Cauchy s 


n(A) = A-f — 
A2 


(A.ll) 


and, X = i!L±i£(M:i 2 T 
, (16n25) 

W ere, s - 1.51 is the refractive index of the 7059 glass and 

mission m the interference free region (fig. A2). 


(A.12) 


Tq is the measured trans- 
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Determination Of Thickness (d): 

In the transparent and medium and weak absorption regions (Fig. A2), we get 
interference fringes. In these regions, d is calculated using the following equation: 

~ 2(Ain2 - Aani) ■ 

where (ni, A^) and (^ 2 , A 2 ) are the values of refractive index and wave length two 
adjacent maxima (or minima), d is calculated for various pair of maxima and minima. 
Average value of = Hdlri) is now calculated. This average value of d is put into 
the equations, 


ndav = mXj2 : for maxima^ 


(A. 14) 


'^dav — (2m + l)A/4 : for minima^ (A.15 i 

where m is an integer, m calculated using Eqs. (14) and (15) is approximately an 
integer. This is made an integer and the new values of n(A) are calculated using Eqs. 
(14) and (15). These new values of n(A ) are put back in Eq. (13) to calculate the 
thickness d and d^v again. This gives the correct value of d. 

New values of n(A) and d are used in Eqs. (10) and (12) to calculate a in the weak 
- medium and strong absorption regions respectively. In the strong absorption region, 
where there is no interference fringes, we extrapolate the calculated n(A) using Cauchy s 
equation given by Eq. (11). 

To determine the optical gap of a-Si:H film {ahi/Y^'^ versus hu is plotted for hu > 
1.5eV. It gives a straight line (see Fig. 2.2). 

{ahiyY^^ = Ax{hy-E,), (A.lbi 

Eg is obtained by the intercept of the straight line. 
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Figure A.l; Schematic showing tramission geometry of absorbing film (a-Si:H) on a 
non absorbing substrate (glass slide) 
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Appendix B 


Movement of the stepper motor 


The stepper motor can be moved clockwise or anti clockwise by applying 12 volts 
across two of the four coils of the motor in a specific sequence. The sequence of the coil 
voltages for clockwise or anti clockwise rotation of the motor is given in Fig. Bl. This 
sequencing is achieved by connecting four electronic switches to the four coils of the 
motor. These switches, when closed, connect a particular lead to the ground voltage 
and hence the current flows through the coil. When the switch is on, the current through 
the coil stops. The motor driver circuit is shown in Fig. B2. The circuit contains four 
opto-isolated relays (OIR, in Fig. B2). An optoisolated relay contains an LED which 
can be switched on by 5volts TTL supplied by the parallel port of the PC. This LED 
in turn switches on a photo transistor and a switch circuit attached with the photo 
transistor. 
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t = 0 At 3At 5At 


0 At 


3At 5At 




CLackwise Rotation 


Anti-clockwise Rotation 


Figure B.l: Voltage sequence for the four coils of the stepper moter for clockwise and 

anticlockwise rotation of the motor. The sequencing is done by the parallel port of the 

PC by sending TTL voltage. At is the delay behveen the two steps and hence controls 
the speed of the motor. 
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12 Volt 



Figure B.2: The stepper motor driv(;r circuit seperate groundings are used for the PC 
and the power supply to safegaurd the PC. Diodes (1N5402) are used to provide low 
resistance path for the quick dissipation of the energy when coil is changed from on to 
off position. 
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Appendix C 

Tbasic program for Transmission and 
PDS 


10 rem programme for PDS/TRANSMISSION using monochromator and motor 
20 rem initialisation of monochromator 
25 els 

27 input ’’print 1 for transmission and 2 for PDS”;r 
30 input ” present wavelength ”;a 
40 input ” starting wavelength ”;b 
45 e = b 
47 t=0 

50 input ’’final wavelength ”;c 
60 input ” step”;d 

65 input ’’enter file. name” ,file.name $ 

67 Open file.name $ for output as #2 

68 print #2. ’’DATE : ”,DATE S,”Time : ”, TIME $ 

70 al = abs (a-b)*5 

75 bl = (abs (b-c)/d) 

78 b2 = abs (660-b)/d 

80 if a>b then goto 1000 else goto 100 

90 rem if b>c then goto 1000 - 
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100 i = 0 
102 rem k = 0 
105 j = 0 

no out khS78, 10 
120 delay .1 
130 i = i+1 

140 if i> al then gosub 2800 

141 rem if t> 0 and cib then goto 1420 

142 if j = b2 then goto 143 else goto 146 

143 input "put filter” ;p 

144 j = j+1 

146 if j>bl then goto 2710 
150 out &:h378, 9 
160 delay .1 
170 i = i+1 

180 if i>al then gosub 2800 

181 rem if t> 0 and cjb then goto 2610 

182 if j = b2 then goto 183 else goto 186 

183 j = j+1 

184 input ’’put filter” ;p 
186 if j>bl then goto 2710 
190 out &;h378, 5 

200 delay .1 
210 i = i+1 

220 if i>al then gosub 2800 

221 rem if t> 0 and cib then goto 1820 

222 if j = b2 then goto 223 else goto 226 

223 j= j+1 

224 input ’’put filter” ;p 
226 if j>bl then goto 2710 
230 out &h378, 6 
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240 delay .1 
250 i = i+1 

260 if i>al then gosub 2800 

261 rem if t 0 and cjb then goto 1820 

262 if j = b2 then goto 263 else goto 266 
263j =j+l 

264 input ’’put filter” ;p 

266 if j>bl then goto 2710 

270 goto no 

1000 i = 0 

1100 out &h378, 10 

1200 delay .1 

1300 i = i+1 

1400 if i> al then gosub 2800 

1420 rem if t>0 and c>b then goto 141 

1462 if j = b2 then goto 1464 else goto 1466 

1466 if j>bl then goto 2710 

1500 out &h378, 6 

1600 delay .1 

1700 i = i+1 

1800 if i>al then gosub 2800 

1820 rem if t> 0 and ob then goto 261 

1862 if j = b2 then goto 1864 else goto 1866 

1864 input ’’put filter” ;p 

1866 if j>bl then goto 2710 

1900 out &h378, 5 2000 delay .1 

2100 i = i+1 

2200 if i>al then gosub 2800 
2220 rem if t>0 and c>b then goto 221 
2262 if j = b2 then goto 2264 else goto 2266 
2264 input ’’put filter” ;p 
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2266 if j>bl then goto 2710 
2300 out &h378, 9 
2400 delay .1 
2500 i = i+1 

2600 if i>al then gosub 2800 

2610 rem if t> 0 and c>b then goto 181 

2620 if j = b2 then goto 2640 else goto 2660 

2640 input ’’put filter” ;p 

2660 if j>bl then goto 2710 

2670 goto 1100 

2680 rem t = t+1 

2710 if r = 2 then goto 2790 

2715 input remove sample and reset the monochromator or print 1 to end”;q 

2720 if q = 1 then goto 2790 

2725 j = 0 

2730 e = b 

2735 rem t = t+1 

2740 goto 70 

2790 end 

2795 goto 5010 

2800 open ”com2:9600,N,8,2,cs,ds,cd” as #1 
2810 print ” ” 

2815 for i = 1 to 200 

2816 next i 

2817 j = j+l 
2828 delay 30 
2820 Print #l,'’q” 

2830 input #l,vl 

2832 if vl > 2e-I then print #I,”g 24”;goto 3000 

2834 if vl< 2^1 and vf> 2e-2 then print #l.”g 23” igoto 3000 

2835 if vKle-l and vl>le-2 then print #l,”g 22”:goto 3000 
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2836 if vl< 5e-2 and vl> 5e-3 then print #l,”g 21”:goto 3000 
2838 if vl<2e-2 and vl> 2e-3 then print #l,”g 20” :goto 3000 
2840 if vl< le-2 and vl> le-3 then print 19” :goto 3000 

2850 if vl<5e-3 and vl> 5e-4 then print #l,”g 18” :goto 3000 
2852 if vl< 2e-3 and vl> 2e-4 then print #l,”g 17”:goto 3000 
2854 if vl<le-3 and vl> le-4 then print #l,”g 16” :goto 3000 
2856 if vl< 5e-4 and vl> 5e-5 then print #l,”g 15”:goto 3000 
2858 if vl<2e-4 and vl> 2e-5 then print #l,”g 14” :goto 3000 
2860 if vl< le-4 and vl> le-5 then print #l,”g 13”:goto 3000 
2862 if vl<5e-5 and vl> 5e-6 then print #l,”g 12”:goto 3000 
2864 if vl< le-4 and vl> le-5 then print #l,”g U”;goto 3000 
2866 if vl<5e-5 and vl> 5e-6 then print #l,”g 10”:goto 3000 
2868 if vl< 2e-5 and vl> 2e-6 then print #l,”g 9”:goto 3000 
2870 if vl<le-5 and vl >le-6 then print #l,”g 8” :goto 3000 
2880 if vl<5e-6 and vl> 5e-7 then print #l,”g 7” :goto 3000 
2890 if vl< 2e-7 then print #1, ”g 6” 

3000 if vl>2e-l then print #l,”tl 7” 

3010 if vl<= 2e-l and vl> 2e-3 then print #1, ”tl 8” 

3020 if vl<= 2e-3 and vl> 2e-6 then print #l,”tl 9” 

3030 if vl<= 2e-6 then print #l,”tl 10” 

3040 delay 30 
3050 print #l,”q” 

3060 input =l,v2:delay 30 
3070 print :^l,”q” 

3080 input #l,v3 
3090 rem print :^l,”q” 

4000 rem input ?5^1,v4 

4005 close #1 

4006 rem k = k-f 
4010 s = (v2+v3)/2 

4020 print ^output at lamda”; s, e 
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4022 rem Open file.name $ for output as #2 
4025 rem Open filename $ for output as #2 
4030 rem Open file.name $ for output as #2 
4045 rem if j = bl then goto 4050 else goto 4060 
4050 rem close ^2 
4060 e = e+d 

4070 rem if j = abs((640-b)/d) then goto 4090 

4080 a = b 

4090 i = 1 :al = d*5 

5000 return 

5010 close #2 

6000 stop 

7000 end 
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Appendix D 


Position sensitive detector (PSD) 


Fig. D1 shows the equivalent circuit diagram of PSD. We use quadrant Si Photodiode 
as PSD. In illuminated condition, a parasitic resistance (Rp) is found between the two 
diode pairs because the diodes are made on the same photosensitive Si wafer. Two 
shunt resistances {Rsk) are used as current to voltage converter as shown in Fig. Dl. 
Voltage signals A and B are obtained through these two diode pairs. 
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Appendix E 


Analog circuit 


Fig. El shows the circuit which amplifies and does the operation (A-B)/(A+B) on 
the outputs A &: B of the PSD. The main part of the circuit consists of two differential 
amplifiers. Output of PSD, A &: B are given in input of one of the differential amphfiers. 
The input of the other differential amphfier is A & (-B). So the signal B is passed through 
an Operational amplifier in the inverter mode as shown in the Fig. A6. Outputs of 
the two differential amplifiers are (A-B) and (A+B) respectively. These signals are 
now passed through Buffer circuits and are fed to a divider IC. Output of the divider 
(A-B)/(A+B) is the PDS signal. Actual components used by us are as indicated in 
Fig. El. 
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Appendix F 


Tbasic program for CPM 

10 els 

12 rem user input parameters 
14 rem 

20 print ’’input parameters” 

30 input ’’file name”, file.name $ 

40 input ’’value of constant photocurrent”;cpc 
50 input ’’present value of lamda”;al 
60 input ’’initial value of lamda”;a2 
70 input ’’final value of lamda”;a3 
80 input ”step”;A4 

90 rem open file.name $ for append as #1 
95 rem j = 0 to abs(a3-a2)/a4 
100 gosub 500 
102 rem print ”ok” 

105 rem end 
107 print ”ok” 

110 for j = 0 to abs(a3-a2)/a4 
120 gosub 1100 

122 print ’’back after reading lock in” 

125 rem end 
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150 gosub 1500 
155 rem. end 
170 gosub 2000 

175 open file.name $ for append as #1 190 close #1 

200 s = a4*5 : i = 0 

210 gosub 502 

215 a2 = a2+a4 

217 rem end 

220 next j 

240 rem close ■ 

250 end 

400 rem subroutine for the motor with monochromator 

401 rem 

500 s = abs(a2-al)*5 :i = 0 

502 if a = 10 then goto 540 

503 if a = 9 then goto 570 

504 if a = 5 then goto 600 

505 if a = 6 then goto 510 
510 out &h378,(10+b) 

515 a = 10 

520 delay .1 : i = i-j-i 
530 if i> s then goto lOOO 
540 out &h378,(9+b) 

545 a = 9 
550 delay .1 ; i = 

560 if i> s then goto 1000 
570 out &h378,(5+b) 

575 a = 5 

580 delay .1 : i = i-j-j 
590 if i > s then goto lOOO 
603 out &h378,(6+b) 
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605 a = 6 

610 delay .1:1 = 1+1 

620 If 1 > s then goto 1000 else goto 510 

1000 return 

1100 rem subroutine for reading the analog lock in amplifier 

1101 rem 

1105 dim dta % (10000) 

1110 port % = &h220 
1120 out port % +11,1 
1130 dummy = inp(port %+4) 

1140 ch % = 0 

1150 rem for Ipl = 0 to 2 

1155 delay 10 

1160 out port % +10, ch % 

1170 for waitl = 0 to 25: next waitl 
1180 out port % +12, 

1190 dh % = inp (port % +5) 

1200 if dh %>15 then goto 1190 
1210 dl % = inp(port % +4) 

1220 dta % = ((dh %*256 +dl % - 2048)''50)/10.23 
1125 print "data = dta % 

1230 rem next Ipl 

1240 rem dta % = (dta %(0) + dta %(1) +dta %(2))/3 
1250 out port % +11,0 
1260 return 

1500 rem comparison with cpm 

1501 rem 

1510 if abs(dta %) i 1.02*cpc and abs(dta %)i.98*cpc then goto 1900 
1520 if abs(dta %) i cpc then goto 1700 else goto 1530 
1530 if b = 160 then goto 1570 
1535 if b = 96 then goto 1600 
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1540 if b = 80 then goto 1630 
1545 if b = 144 then goto 1550 
1550 out &h378, (160+a) :b = 160 
1560 delay 1 : gosub 1100 
1562 rem print ’’hello” 

1565 goto 1510 

1567 rem print ’’done 

1570 out &:h378, (96+a) :b = 96 

1580 delay 1 : gosub 1100 

1585 goto 1510 

1600 out &h378, (80+a) ;b = 80 
1610 delay 1 : gosub 1100 
1620 goto 1510 

1630 out &:h378, (144+a) ;b = 144 
1640 delay 1 : gosub 1100 
1650 goto 1510 

1700 rem if dta %> cpc then goto 1530 else goto 1700 

1710 if b = 160 then goto 1760 

1715 if b = 144 then goto 1790 

1720 if b = 80 then goto 1820 

1725 if b = 96 then goto 1730 

1730 out &;h378, (160+a) :b = 160 

1740 delay 1 : gosub 1100 

1750 goto 1510 

1760 out &h378, (144+a) ;b = 144 
1770 delay 1 : gosub 1100 
1780 goto 1510 

1790 out &h378, (80+a) ;b = 80 
1800 delay 1 : gosub 1100 
1810 goto 1510 

1820 out &h378, (96+a) ;t = 96 
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1830 delay 1 : gosub 1100 
1850 goto 1510 
1900 return 

2000 rem subroutine to read SR510 

2001 rem 

2010 open ”com2:9600,n,8,2,cs,ds,cd” as #2 
2020 print #2,” ” 

2030 for i = 1 to 200: next i 
2040 print #2, ”br’:delay 1 
2050 print #2,”q” 

2060 input #2,vl 


2065 

rem print #2,”( 

1” : 

dela} 

1 . 





2070 

if 

vl< 

2e-l 

and 

vl>2e-2 

then 

print 

#2,” 

g23” 

:goto 2270 

2080 

if 

vl< 

le-1 

and 

vl> 

le-2 

then 

print 

#2,’ 

>g22’ 

-•goto 2270 

2090 

if 

vl< 

5e-2 

and 

vl> 

5e-3 

then 

print 

#2,’ 

’g2r 

:goto 2270 

2100 

if 

vl< 

2e-2 

and 

vl> 

2e-3 

then 

print 

m: 

’g20’ 

:goto 2270 

2110 

if 

vl< 

le-2 

and 

vl> 

le-3 

then 

print 

m: 

’g 19’ 

:goto 2270 

2120 

if 

vl< 

5e-3 

and 

vl> 

5e-4 

then 

print 

m: 

’g 18’ 

:goto 2270 

2130 

if 

vl< 

2e-3 

and 

vl> 

2e-4 

then 

print 

#2,’ 

’g 17’ 

:goto 2270 

2140 

if 

vl< 

le-3 

and 

vl>le-4 1 

then print : 

#2,” 

gl6” 

:goto 2270 

2150 

if 

vl< 

5e-4 

and 

vl> 

5e-5 

then 

print 

m: 

’g 15’ 

:goto 2270 

2160 

if 

vl< 

2e-4 

and 

vl> 

2e-5 

then 

print 


>g i4>^ 

:goto 2270 

2170 

if 

vl< 

le-4 

and 

vl> 

le-5 

then 

print 

m: 

’g 13’^ 

:goto 2270 

2180 

if 

vl< 

5e-5 

and 

vl> 

5^6 

then 

print 

m: 

’g 12’^ 

:goto 2270 

2190 

if 

vl< 

2e-5 

and 

vl> 

2e-6 

then 

print 

m: 

’g ir 

:goto 2270 

2200 

if 

vl< 

le-5 

and 

vl> 

le-6 

then 

print 

#2,’' 

’g 10’^ 

:goto 2270 

2210 

if 

vl< 

5^6 

and 

vl> 

5e-7 

then 

print 

#2,’' 

’g 9” 

:goto 2270 

2220 

if 

vl< 

2e-6 

and 

vl> 

2e - 1 

then 

print 

#2/ 

’g 8” 

.•goto 2270 

2230 

if 

vl< 

le-6 

and 

vl> 

le-7 

then 

print 

#2,’ 

>g 7” 

:goto 2270 

2240 

if 

vl< 

5e-7 

and 

vl> 

5e-8 

then 

print 

#2,’ 

>g 6” 

:goto 2270 

2250 

if 

vl< 

2e-7 

and 

vl> 

2e-8 

then 

print 

#2/ 

-g 5” 

:goto 2270 
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2260 rem if vl< 2e-l and vl> 2e-2 then print #2,”g 23” 
2270 rem if vl< 2e-l and vl> 2e-2 then print #2,”g 23” 
2280 if vl>= le-2 then print #2,”tl, 8” 

2290 if vl< le-2 and vl>= le-3 then print ^2,”tl, 8” 
2300 if vl< le-3 and vl>= le-5 then print #2,”tl, 9” 
2310 if vl< le-5 then print ^2,”tl, 9” 

2320 delay 15 
2330 print 5i^2,”q” 

2340 input #2, v2 
2350 delay 5 
2360 print #2,”q” 

2370 input #2, v3 
2380 delay 5 
2390 d = (v2 -1- v3)/2 
2400 rem input #l,d 
3000 close 
3100 e = a2 

3200 print ”output at lamda ”:d, e 
4000 return 
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